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(54) Title: METHOD FOR CONTROLLING CHEMICAL SYNTHESIS PROCESSES 
(54) litre: METHODE DE REGULATION DE PROCEDES DE SYNTHESE DE PRODUTTS CHMQUES 
(57) Abstract 

A method for controlling a process for synthesis ing at least one chemical in 
a plant including at least one reactor (R) mat may be considered as a perfectly 
mixed reactor, wherein control variables (GC) enable the process to be controlled 
so that a number of variables known as controlled variables (GR) and related to the 
properties of the product and/or the operation of the process match corresponding 
set values (Cor). Said method comprises the steps of (a) inputting set values for the 
controlled variables (Cqr); (b) using a prediction unit (OP) to calculate predictions 
of the controlled variables (Pan) on the basis of measurements of die process control 
variables (Mqc); (c) using a control unit (OC) to calculate set values for the process control variables (Cqc) on the basis of die set values 
(Cqr) and predictions (Pgr) of the controlled variables; and (d) transmitting set values for die process control variables (Coc) to process 
control actuators. The prediction unit (OP) is based on a mathematical process model known as a direct model (MX and designed in such 
a way mat the mass Mxr of at least one component (X) in the reactor (R) is predicted by means of die equation Mxr - LAG (FxRfat * r% 
rx), wherein FxRia is the mass flow rate of component X at the inlet of the reactor R, rx is the dwell time of x in die reactor, and the 
function y - LAG (u, r) is the solution of die differential equation r . dy/dt + y, calculated with die instantaneous value of u and r as 
well as with die last calculated value of y. 

(57) Abrege* 

Une mCthode de regulation d'un precede de synthese d'au moms un produit chimique dans un appareillage comprenant au moins un 
reacteur (R) pouvant etre assimile* a un reacteur parfaitement melange, dans laquelle des grandeurs de commande (GC) permettent d'agir sur 
le deroulement du proceed en vue que des ou plusieurs grandeurs liees aux proprietes du produit et/ou au deroulement du precede, appelees 
grandeurs r£ glees (GR), soient egales a des consignes correspondantes (Cqr), ladite m 6th ode comprenant les etapes suivantes: (a) entree 
de consignes concemant les grandeurs rfglees (Cqr); (b) calcul, au moyen d'un organe de prediction (OP), de predictions des grandeurs 
nSglees (Pgr), sur la base de mesures des grandeurs de commande du urocele* (Moc); (c) utilisation d'un organe de contrdle (OC) pour 
calculer des consignes des grandeurs de commande du proceed* (Cgc)» sur base des consignes (Cqr) et des predictions (Pgr) des grandeurs 
reglees; (d) transmission des consignes des grandeurs de commande du precede (Cgc) & des actionneurs afin d'agir sur le deroulement du 
proo&te; dans laquelle 1' organe de prediction (OP) est base sur un modele mathematique du precede, appeie modele direct (M),*t est concu 
de telle maniere qu'on pnSdit la masse Mxr d'au moins un constituent (X) dans le reacteur (R) par l'equation: Mxr ■ LAG (FxRin . rx. 
rx), dans laquelle FxRin est le debit massique du constituant X entrant dans le reacteur R; rx est le temps de sejour de x dans le reacteur, 
la fonction y - LAG (u, r) est la solution de l'equation differentielle r . dy/dt + y t calculee avec la valeur instantanee de u et de f, atnsl 
qu'avec la demiere valeur de y calculece. 
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M&hode de regulation de precedes de svnthese de produits chimiques 



La prdsente invention concerne une methode de regulation de 
procedes de synthese de produits chimiques. EUe concerne egalement un 
dispositif de regulation pour la mise en oeuvre de cette methode, ainsi qu'un 
procede de synthase, en particulier de polymfere, regule par cette methode. 
5 Dans un precede de synthase de produits chimiques conduit 

conventionnellement, des regulators de type PID (proportionnel - integral - 
diffftrentiel) sont utilises pour reguier individuellement un nombre plus ou 
moins important de grandeurs (temperatures, debits, pressions, ...) ayant 
une influence sur le deroulement de la synthase. En d'autres termes, pour 

10 chaque temperature, debit ou pression a reguier, on mesure en continu (ou 
par intermittence) sa valeur effective, et un regulateur PID compare cette 
valeur effective h une valeur de consigne et agit sur la grandeur a 
commander de fa$on a reduire, le cas echeant, la difference entre la valeur 
de consigne et la valeur mesurge. 

15 Vu la complexity de la plupart des precedes industriels de synthase 

chimique, les valeurs de consignes des differents regulateurs doivent 
aujourd'hui encore fetre ajustees empiriquement pour obtenir fmalement les 
propriety souhaitees du produit synthdtise. On utilise a cette fin des 
recettes, qui fournissent des combinaisons de parametres determines 

20 empiriquement pour obtenir, en regime stable, les proprietes souhait6es du 
produit synthetise. 

De ces recettes pounaient etre deduites, k 1'aide d'outils statistiques 
plus ou moins sophistiques, des relations empiriques entre les grandeurs 
regimes et les proprietes du produit synthetise. H est cependant evident que 

25 ces relations empiriques ne peuvent guere tenir compte des multiples 

interdependances existant entre les differentes grandeurs r6glees separement, 
pas plus que des perturbations inconnues telles que les teneurs en impuretes 
des matieres premises. 

n est aussi Evident qu'une regulation classique k boucle fermee, uti- 

30 lisant comme corrections de retour des mesures de proprietes essentielles du 
produit synthetise, est difficilement applicable a la plupart des procedes de 
synthese. En effet, les temps morts intervenant soit dans le procede, soit 



dans les mesures ou analyses utilisees comme corrections de retour, sont 
trop eleves, et les intendependances entre les differentes grandeurs regissant 
le precede sont trop complexes. 

II est connu depuis longtemps que les procedes de synthese, 
notamment les procedes de synthese en continu de polymeres (procedes de 
polymerisation), equipes de regulations avec ajustage empirique des valeurs 
de consignes, ont des inconvenients majeure, qui peuvent etre resumes 
, comme suit : 

- les demanages du precede de synthese prennent beaucoup de temps et 
geherent des quantity importantes de produit "hore normes" ; 

- les transitions de qualite sont lentes, ce qui entrame egalement la produc- 
tion de quantites importantes de produits de transition "hore normes" ; 

- I 'allure de marche du precede, e'est-a-dire la masse de produit(s) 
syntfaetisee par unite de temps, peut difficilement etre modifiee sans 
alterer les propridtes de ce ou ces produits ; 

- la constance des proprietes essentieUes du ou des produits synthases 
laisse souvent a desirer, meme en regime stable. 

Afin d'eviter un ajustage empirique des valeurs de consigne, il a ete 
propose dans la literature specialisee d'utiliser des methodes de regulation 
de proceeds de synthese qui ont recoure a des equations caractdristiques 
moddlisant le proeddd de synthese pour.relier certaines propridtds du ou des 
produits synthdtisds aux conditions de marche du ou des reacteure lore de la 
synthese. Cependant, afin de limiter la complexity de ces equations caracte- 
ristiques, on considerait jusqu'a present qu'on dtait oblige en pratique soit 
de considdrer exclusivement le cas statique (regime stable), soit de se 
limiter a une modelisation empirique tres simplifiee de la dynamique du 
proedde. L'utilisation d'un modele statique est limitee au controle d'un 
regime de production assez stable. 

Dans le cas d'une modelisation empirique, les equations 
caractdristiques ne sont valables que pour une plage de validity etroite (a 
proximhd du point ou la modetisation a dtd effectuee). Dans les deux cas, 
les phases de demarrage et les phases de transition sont mal maitrisees. II 
est certes envisageable de "couvrir" une plage plus large de conditions de 
fonctionnement en procedant a une telle modelisation locale en plusieure 
points distincts de 1'espace des parametres de fonctionnement, mais une telle 
approche devient prohibitive des que Ton tente de reguler plusieure gran- 
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deurs en agissant sur de multiples param&res. 

II serait par consequent souhaitable de disposer d'une methode et 
d'un dispositif de regulation simples, mieux adapts aux specifiers de la 
dynamique de precedes de synthase de produits chimiques. 

A cette fin, ia pi^sente invention concerne une methode de regulation 
d'un precede de synthese d'au moins un produit chimique dans un 
appareillage comprenant au moins un reacteur (R) pouvant etre assimiie k 
un reacteur parfaitement melange, dans laquelle une ou plusieurs grandeurs 
de commande (GC) permettent d'agir sur le deroulement du proc&te en vue 
qu'une ou plusieurs grandeurs li&s aux proprietes du produit et/ou au 
deroulement du precede, appeiees grandeurs reglees (GR), soient egales a 
des consignes correspondantes (C QR ) (ou du moins aussi proches que 
possible de celles-ci), ladite metbode comprenant les etapes suivantes : 
(a) entrfe de consignes conceniant les grandeurs reglees (C GR ); 
15 (b) calcul, au moyen d'un organe de prediction (OP), de predictions des 
grandeurs reglees (Pq R ), sur la base de mesures des grandeurs de 
commande du precede (M GC ); 

(c) utilisation d'un organe de contrdle (OC) pour calculer des consignes 
des grandeurs de commande du precede (C GG ), sur base des 

20 consignes (C GR ) et des predictions (Pq R ) des grandeurs reglees; 

(d) transmission des consignes des grandeurs de commande du precede 
(Cqc) * des actionneurs, ou a des organes de reglage contrdlant des 
actionneurs, afm d'agir sur le deroulement du precede; 

dans laquelle l'organe de prediction (OP) est base sur un modfele 
25 matheniatique du precede, appeie modfele direct (M), et est con$u de telle 
maniere qu'on predit la masse M^r d'au moins un constituant (X) dans le 
reacteur (R) par requation : 

= LAG (Fjq^ . r x , t x ) 

dans laquelle : 

- F XRin est le massique du constituant X entrant dans le reacteur R ; 

- r x est le temps de sejour de x dans le r&tcteur (constante dei temps), qui 
vaut 

7 X = M XR 1 & Fxdis > 

dans lequel : 

M XR designe la derniere valeur calculee de la masse du 
constituant X presente dans le reacteur R ; 



30 
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- E Fxdis d&igne la somme de tous les debits riiassiques Fxdis avec 
iesquels le constituant X disparait du reacteur R, notamment par 
reaction et/ou par sortie du reacteur ; 
- la fonction y = LAG (u , r) est la solution de l'equation differentieUe 

u = t • + y 

calculee avec la valeur instantanee de u et de r, ainsi qu'avec la derniere 
valeur de y calcutee. 

L'avantage de cette m&hode est que Tequation differentielle 
susmentionnee est resolue par un calcul algebrique simple, par exemple par 
la formule suivante (T designant Pintervalle de temps, generalement petit 
par rapport a r, separant les calculs successifs) ou par une fonnule 
equivalente a celle-ci : 

y(t-T) * u(t)-— I— 

y(t) = ; *(t) 

i ♦ _L_ 

T (t) 

Dans le cas ou les masses de plusieurs constituants sont evaluees 
comme expose ci-dessus, la m&hode de rinvention est paiticulierement 

avantageuse dans la mesure oil ces masses peuvent 6tre calculees 

15 sequentiellement par de tels calculs algebriques simples, recalcules 
frequemment (T < < r en g&uSral). Par opposition, les methodes 
traditionneUes necessitent la resolution simultanee d'un systeme d'equations 
differentielles, ce qui exige gdn^ralement une grande puissance de calcul et 
des algorithmes de calcul numenque (integration) sophistiqufc ; il en resulte 

20 que la duree de chaque iteration de calcul est elevee, et par consequent une 
regulation de ce type reagit mal a des variations rapides. 

Le proce<ie de synthese regule peut servir a la synthase d'un 
compost monomerique ou polymerique ; de ties bons resultats ont 6t6 
obtenus pour la regulation de precede* de polymerisation. Le precede 

25 s'etend aussi au cas ou plusieurs produits utiles sont synthetises 

simultanement dans le meme preceded . Le procette peut etre continu ou 
discontinu (batch) ; la methode de regulation de rinvention donne 
d'excellents results dans le cas de precedes continus. Le precede de 
synthese regule peut eventueUement ne constituer qu'une partie d'un 

30 procede plus vaste, dont les autres parties sont regulees d'une maniere 

differente, ou ne sont pas regulees. Pour que la methode de regulation de 



l'invention soit applicable, il faut qifau moins un reacteur puisse etre 
assimild a un reacteur parfaitement melangd, c'est-a-dire a un i^acteur dans 
lequel les difterentes grandeurs (temperature, concentrations des constituants 
presents, etc.) sont quasiment identiques en chaque point. D'autres 
i^acteurs eventuels peuvent etre du type piston ("plug-flow") : on les 
modelise mathematiquement par des temps morts. La methode est 
dgalement applicable a un procedS se ddroulant dans plusieurs rtacteurs 
disposes en serie et/ou en paralfele, pouvant produire des produits aux 
propri&6s identiques ou differentes. 

Par "constituants", on entend designer rensemble des substances 
prtsentes dans le reacteur et destindes a paiticiper St la synthase ou k la 
permettre : non seulement les reactifs de depart ainsi que le ou les produits 
synth£tis6s, mais €galement la ou les eventuelles substances ne subissant 
aucune transformation, telles que solvants, catalyseurs, etc. 

Outre un ou plusieurs reacteurs, V installation dans iaquelle se 
d&oule le proc&te peut eventuellement comprendre d'autres dispositifs 
classiques tels que d&endeurs, strippers, condenseurs, secheurs, colonnes de 
distillation, etc. Generalement, ces dispositifs auxiliaires peuvent d'ailleurs 
etre egalement consid6r6s comme des reacteurs (parfaitement melanges ou 
de type piston), meme si aucune reaction chimique ne s'y deroule. 

Dans le cas d'un procddd de polymerisation, les "grandeurs liees aux 
proprietes du produit" peuvent par exemple 6tre choisies panni la masse 
mol&ulaire, Findice de fluidite en fondu, la masse volumique standard, la 
teneur en comonomere lorsqu'un comonomere est prtsent, etc. 

Des exemples de "grandeurs lides au procedS" sont notamment la 
temperature et la pression regnant dans le reacteur, Failure de marche du 
proc&te, les concentrations des differents reactifs dans le reacteur, etc. 
L' allure de marche designe la masse de produit synth6tis6e par unit£ de 
temps, qui n'est cependant pas forc&nent egale au ddbit de produit 
synth&ise sortant du reacteur : ainsi, & titre d'exemple, dans les phases de 
demarrage notamment, le debit massique de produit synth&ise sortant du 
reacteur est generalement tres faible, voire nul, bien que la synthese ait 
commence, c'est-a-dire que ce debit sortant est alors inf&ieur a Failure de 
marche. En regime stable, par contre, on peut assimiler Failure de marche 
a la masse de produit synthetise par unit6 de temps. 

Des exemples de "grandeurs de commande" sont les debits de 
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reactifs entrant dans le reacteur, la puissance foumie aux dispositifs de 
chauffage, etc. Ce sont ces grandeurs qui pemiettent d'agir sur le 
deroulement du precede ainsi que, generalement, sur les proprietes du 
produit synthetise. 

5 La ou les consignes de la ou des grandeurs de commande sont 

transmises directement ou indirectement a des actionneun! classiques tels 
que notamment des vannes, des elements chauffants, etc. "Indirectement" 
signifie que les grandeurs de commande peuvent etre transmises par 
1'intennediaire d'un ou plusieurs organes de regulation (controlant 
10 generalement une seule variable, par ex. des regulateurs HD) controlant le 
ou les actionneurs (regulation "locale"). 

Sur le plan materiel, l'organe de prediction et 1'oigane de controle 
sont generalement des dispositifs de calcul classiques permettant d'effectuer 
des calculs en fonction de leur cablage ou de leur programmation ; il peut 
15 notamment s'agir d'ordinateurs ou de systemes numeriques de controle- 

commande (SNCC). Un seul dispositif peut avantageusement combiner les 
fonctions de prediction et de controle. Le ou les dispositifs de calcul 
utilises sont de preference du type numerique, et fournissent pdriodiquement 
(par intermittence) les resultats de leurs calculs. Les intervaUes de temps 
20 separant la foumiture de ces resultats peuvent varier dans le temps, et 

peuvent egalement differer selon le resultat conceme : il est clair que des 
grandeurs a variation rapide doivent etre recalculees plus frequemment que 
des grandeurs k variation lente. Des registres a decalage peuvent etre 
utilises pour simuler materiellement les temps morts. 
25 L'organe de prediction est base sur un modele math&natique direct 

du pwc6d6 (M), dans lequel le reacteur (R) est assimile k un reacteur 
parfaitement melange ; un ou plusieurs d<Slais purs (temps morts) peuvent 
eventuellement etre pris en compte pour representer d'6venmels reacteurs de 
type piston, d'eventuels delais de transport ou delais d'obtention de r&ultats 
30 de mesures, .... 

L'organe de controle est de preference base sur Finverse du modele 
direct utilise dans l'organe de prediction (modele inverse). 

En general, la somme E Fxdis de tous les debits massiques (Fxdis) 
avec lesquels le constituant X disparait du reacteur R comprend deux 
35 termes : 

* FR X qui designe le debit massique suivant lequel X est consomme par 
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une ou plusieurs reactions chimiques eventuelles ; 
* F Xout 9 ui d &igne reventuel debit massique de X sonant du reacteur 
par soutirage en cours de reaction, dans le cas (habituel) oil X n'est pas 
entifcrement consomme par reaction dans ce reacteur ; ou encore, par 
5 exemple, par evaporation dans le cas d'un reacteur ouvert. 

L'interet de la methode est que les termes Fxdis sont en general 
proportionnels a ; par exemple, on a en general 
F Xout = M XR L r R 

(r R ddsignant le temps de s6jour du reacteur R) 
10 et Frx = R x . MxR 

(R x designant la reactivite de X dans le reacteur R). 

Dans ce cas, Texpression donnant r x se simplifie et devient : 

r x = I / (R x + 1/ TR ) 

Cette expression est independante de M^, ce qui constitue une 
15 simplification extremement interessante. 

Un autre avantage de la methode reside dans le calcul periodique du 

temps de sejour r x . En effet, t x repr&ente bien la dynamique du 

constituant considers dans le reacteur. Ceci permet notamraent de suivre 

Involution de ce parametre, ce qui est important pour la comprehension de 
20 la dynamique du precede, et par consequent pour sa regulation. Au 

contraire, les methodes empiriques de type "boite noire" ne permettent pas 

d'acceder a ce parametre. 

Avantageusement, le calcul des predictions des grandeurs regimes 

(PqR^ P^ 1 additionnellement tenir compte d'une ou plusieurs mesures de 
25 grandeurs r6g!6es (Mq R ), de grandeurs de commande QAqq) et/ou d'autres 

grandeurs liees au deroulement du precede 

De meme, avantageusement, le calcul des consignes des grandeurs 

de commande du precede (C GC ) peut additionnellement tenir compte d'une 

ou plusieurs mesures de grandeurs iegiees (M GR ), de grandeurs de 
30 commande (M GC ) et/ou d'autres grandeurs liees au deroulement du precede 

^AP)' iden *iques ou differentes de celles eventuellement prises en compte 

pour le calcul des predictions des grandeurs rfgiees (Pq R ). 

Toutes les mesures dont il est question dans la presente description 

ne sont pas forcement des mesures directes, dans le sens oft Tune ou 
35 plusieurs d'entre elles peuvent eventuellement etre des mesures 

inferentielles, e'est-a-dire des valeurs obtenues par calcul a partir d'une ou 
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plusieurs autres mesures directes. Ainsi, par exemple, Tallure de marche 
de certains precede* de synthese exothermiques ne peut pas se mesurer 
directement, mais on peut en obtenir une mesure inferentielle par calcul, par 
exemple a partir de mesures (directes) du deoit et des temperatures d'entree 
5 et de sortie du fluide de refroidissement. 

Dans le cas particulier de precedes de polymerisation, la ou les 
propriety du polymere intervenant dans la regulation sont de preference 
choisies parmi la masse volumique specifique (MVS) du polymere, les 
proprietes rheologiques du polymere a l'etat fondu, et sa teneur en 
10 comonomere. En particulier, la ou les proprietes rheologiques intervenant 
dans la methode de regulation sont avantageusement 1'indice de fluidity en 
fondu (melt index) du polymere et/ou une mesure de viscosite. 

Ayantageusement, on evalue une ou plusieurs proprietes du polymere 
en utilisant une technique choisie parmi la spectroscopic infrarouge proche 
15 (MR), la spectroscopic infrarouge a transformee de Fourier (FUR) et la 
resonance magnetique nucleaire (NMR). 

En particulier, on peut avantageusement evaluer une ou plusieurs 
proprietes du polymere en appliquant une relation de correlation pre^tabUe 
aux resultats de mesures reaUsees par spectroscopic infrarouge proche (NIR) 
a plusieurs longueurs d'onde pretletenninees en fonction de la nature du 
polymfere, choisies entre 0,8 et 2,6 mm. 

Davantage de detaUs concernant la realisation de telles mesures dans 
le cadre de la regulation de process de polymerisation peuvent etre trouves 
dans la demande de brevet EP 328826 (US 5155184). 
25 En vue de tenir compte d'ecarts eventuels entre les mesures et les 

predictions des grandeurs reglees, il peut etre utile de recourir a une 
collection. 

Un premier type de correction consiste en ce que la consigne d'au 
moins une grandeur reglee (C GR ) soit corrigee sur la base de l'ecart 
(avantageusement filtr6) entre la mesure (M GR ) et la prediction (P GR ) de 
cette grandeur reglee, de facon a ce que la regulation soit efficace (M QR = 
C GR ) meme en presence d'une erreur dans la prediction de cette grandeur 
reglee. Cette technique est couramment designee par l'expression "internal 
model control" (IMC). 

Un second type de correction consiste en ce que le modele (M) du 
prrcede soit periodiquement adapte sur la base de l'ecart (avantageusement 
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filtre) entre les predictions (Pqr) et les mesures (M QR ) des grandeurs 
rggiees, de fa^on a ce que, ici egalement, le modele du precede fournisse 
des predictions des grandeurs reglees (Pqr) aussi proches que possible 
(id&tlement, egales) des mesures de ces grandeurs (Mq R ), ce qui est 
5 indispensable a une regulation efficace. 

L'adaptation consiste a recalibrer le modfcle, c'est-k-dire a en 
recalculer un ou plusieurs param&res ; normalement, le nombre de 
■ parametres recalcutes ne depasse pas le nombre des grandeurs reglees pour 
lesquelles on dispose k la fois d'une prediction et d'une mesure. Une 

10 ^synchronisation (recalage dans le temps) de ces mesures est sou vent 
souhaitable, sunout lorsqu'il s'agit de mesures de proprietes du produit 
synth^tisd dont la duree detention est longue. Ce second type de 
correction est plus avantageux dans la mesure bu il permet d'adapter le 
modele egalement sur le plan de sa dynamique. 

15 L'adaptation concerne non seulement le modele direct du precede 

(organe de prediction), mais egalement le modele inverse (organe de 
controle). 

Selon une variante avantageuse, les mesures (Mq R ) des grandeurs 
r£gl£es interviennent uniquement dans l'adaptation eventuelle du module du 
20 precede, et n' interviennent pas directement dans le calcul des consignes des 
grandeurs de commande du precede (Cq C ). 

C'est-a-dire que les mesures des grandeurs reglees n' interviennent 
pas dans la regulation proprement dite : l'avantage en est que la qualite de 
la regulation n'est ainsi pas affectee par la lenteur eventuelle de revaluation 
25 des proprietes du produit. 

Un autre aspect de l'invention concerne une methode de regulation 
telle que d&rite ci-dessus, appliquee & un precede de polymerisation, 
comprenant Tune ou plusieurs des etapes suppldmentaires suivantes : 
- calcul d'une consigne de temperature dans le reacteur en fonction d'une 
30 ou plusieurs consignes de proprietes du produit ; et transmission de 

cette consigne de temperature a un ou plusieurs actionneurs permettant 
de modifier la temperature dans le reacteur (eventuellement - 
indirectement, c'est-a-dire par rintermediaire d'un ou plusieurs organes 
de regulation, par ex. des regulateurs PID, controlant le ou les 
35 actionneurs) ; 

calcul d'un bilan thermique pour le reacteur, sur base notamment de 
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mesures de temperature; utilisation de ce bilan thermique de fa$on k 
determiner la quantite de polymfere synthetise par unite de temps (allure 
de marche) et/ou la productive du catalyseur et/ou la concentration 
d'au moins un reactif dans le reacteur ; 
5 - calcul de la quantity de chaleur produite par la polymerisation, grace a 
un calcul de la quantite du ou des reactifs qui polymerisent ; 
determination par ce biais de la quantite la chaleur qu'il faut ajouter ou 
evacuer pour maintenir la temperature du reacteur ; utilisation du 
resultat dudit calcul (par exemple par feed-forward) pour ameiiorer la 
10 regulation de temperature, de fagon k respecter au mieux la consigne de 

temperature, notamment en cas de modifications de Failure de marche. 
Ces variantes sont basees sur le lien qui existe entre la quantite du 
ou des reactifs participant k la reaction et la quantite de chaleur produite ou 
absorbee par la reaction. 
15 Selon une variante avantageuse, la propriete Px R d'un constituant 

"x" dans le reacteur R, assimiie a un reacteur parfaitement melange, est 
calculee comme suit : 

Px R = LAG (Px m , Mx R / Fx m ) 
ob "Px" est une propriete d'un constituant M x\ repondant de maniere 
20 substantielle k la loi de melange lineaire Pxj +2 = wj . Pxj + w 2 . 

Px 2 , 

W| et w 2 etant Ies proportions massiques de deux fractions 1 et 2 de 

propriete Px| et Px 2 que Ton melange (wj + w 2 = 1) ; 
Pxj +2 est la propriete de x k sa sortie du reacteur, aprfes melange ; 
25 ^EN est la P ro Priete du constituant M x H a son entiee dans le reacteur R ; 
Mxjj est la masse du constituant x dans le reacteur R ; 
Fx IN est le d&)it m assique du constituant x entrant dans le reacteur R. 

Une transformation mathematique permet parfois de rendre lineaires 
(additives) certaines grandeurs.qui ne le sont pas : par exemple, Tindice de 
30 fluidite en fondu d'un polymere (melt index) ne nSpond pas a une loi de 
melange lineaire, mais bien son logarithme ; le calcul de Px| + 2 
susmentionne s'effecue des lors sur le logarithme de ce parametre. 

Selon une autre variante avantageuse, la methode de regulation de 
T invention comprend les etapes suivantes : 
35 - entree de consignes relatives k une ou plusieurs proprietes du produit a 
• synthetiser, dans un algorithme maitre ; 



- entrfe de la consigne d 'allure de marche du proc6d6; dans un algorithme 
esclave; 

- calcul des consignes de concentration des constituants dans le rgacteur k 
l'aide de 1'algorithme maitre, en fonction notamment des consignes et des 
mesures de propriety du produit ainsi que de mesures ou de predictions 
des concentrations des differents constituants dans le rdacteur ; 

- transmission des consignes de concentration calcul&s par 1'algorithme 
maitre comme grandeurs d'entrfe k 1'algorithme esclave; 

- calcul de consignes de debit des constituants entrant dans le i^acteur, k 
Taide de 1'algorithme esclave, en fonction notamment de la consigne 

d 'allure de marche du precede, des consignes de concentration, et de 
mesures de debit des constituants entrant dans le reacteur, et 

- transmission des consignes de debit calculees a l'aide de 1'algorithme 
esclave k un ou plusieurs actionneurs (eventuellement indirectement, 
c'est-i-dire par 1'intermediaire d'un ou plusieurs organes de rtgulation, 
par ex. des regulateurs PID, controlant le ou les actionneurs) afin de 
regler les debits des constituants entrant dans le reacteur, 

dans Iaquelle 1'algorithme maitre et/ou 1'algorithme esclave sont mis en 
oeuvre comme d&rit ci-dessus, c'est-*-dire en utilisant la fonction LAG 
pour calculer la masse d'au moins un constituant dans le reacteur. 

Les algorithmes maitre et esclave sont egalement mis en oeuvre au 
moyen d'un ou plusieurs dispositifs de calcul classiques. Selon une variante 
avantageuse, I'ensemble des calculs (prediction, contrile, etc.) de ces deux 
algorithmes sont effectuds par un meme dispositif de calcul. 

Avantageusement, des mesures de temperature (par ex. temperature 
dans le reacteur, et/ou temperature d'entree et/ou de sortie d'un eventuel 
fluide de refroidissement) interviennent comme grandeurs d'enti^e 
supplementaires dans l'organe de prediction et/ou de contrdle. 

De preference, 1'algorithme esclave prend egalement en compte des 
mesures de la composition de constituants presents dans ou sortant du 
reacteur. 

De maniere avantageuse, la m&hode de regulation comprend en 
outre une etape de calcul k l'aide de 1'algorithme esclave, en fonction des 
mesures de debits, de predictions de concentrations transmises a 
1'algorithme maitre pour calculer des predictions de proprietes servant 
comme grandeurs d'entree supplementaires dans le calcul des consignes de 
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concentrations. 

Les algorithmes maitre et esclave ferment une regulation du type 
cascade. II est particulierement avantageux que I'algorithme maitre et/ou 
I'algorithme esclave soient adaptatifs, c'est-a-dire que certains de leurs 
5 parametres soient penodiquement recalcules (a intervalles reguiieni ou 
irregubers). Une telle adaptation permet notamment de garantir que le 
modele mathematique decrit le plus fidelement possible le procAJe dans son 
etat courant, meme en cas de modification de certaines conditions 
operatoires (temperature, pression, aUure de marche, etc.) et en cas de 
10 perturbations (empoisonnement du catalyseur, ...). 

L'algorithme maitre realise le controle des proprietes du produit a 
1'aide d'un modele base sur des equations caracteristiques reuant les 
proprietes du produit aux concentrations des difrerenrs constituants dans le 
reacteur, ainsi qu'eventueUement a la temperature regnant dans le reacteur. 
15 L'algorithme esclave regule les concentrations d'un ou plusieurs constituants 
en agissant sur les debits d'alimentation d'un ou plusieurs constituants 
eventuellement differents. 

L'interSt de cette cascade "maitre-esclave" reside dans le fait que le 
modele maitre determine avec precision les concentrations de constituants 
20 necessaires a l'obtention des proprietes voulues pour le produit synthetise, et 
que le modele esclave assure le respect des valeure imposes par le maitre. 
Pilots par le maitre, 1 'esclave est par consequent capable : " 
- d'amener rapidement les concentrations aux valeurs d&irees par le maitre 
et de les y maintenir; 
25 . - de controler efficacement 1'aUure de marche du precede sans perturber 
les concentrations. 

Cette cascade maitre-esclave est particulierement efficace car aussi 
bien le maitre que l'esclave prennent en compte la dynamique du precede 
grace a I'utilisation de la fonction LAG dans les calculs. 
30 L'algorithme esclave peut en outre etre concu de maniere a fournir a 

1'algoritlune maitre des predictions de concentrations fiables. De ces 
predictions ou de mesures de concentrations, ralgorithme maitre deduit des 
predictions de proprietes fiables du produit en train d'etre synthetise dans le 
reacteur. En comparant ces predictions de propriety aux consignes de 
35 proprietes, ralgorithme maitre peut. le cas echeant, intervenir et corriger 
les consignes de concentrations. Cette correction est possible avanr meme 
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qu'un ecart entre une variable et sa consigne ne se presente. La prise en 
consideration de predictions de proprietes obtenues & partir de predictions 
ou de mesures de concentrations pennet de reduire considerablement les 
fluctuations temporelles des propriety du produit synthetise et il en resulte 
. 5 une meilleure Constance de la quality du produit. 

Si les propriety du produit & synth&iser sont d^endantes de la 
temperature dans le ou les r£acteurs, il est preferable de pievoir une 
regulation de la temperature par l'algorithme esclave. Ce dernier dresse le 
bilan thermique de chaque reacteur et determine, au moyen du calcul de 
10 r allure de marche, la quantite de chaleur qu'il faut ajouter ou evacuer pour 
respecter des consignes de temperature calcuiees par l'algorithme maitre. 
De ces resultats il deduit des consignes d'entree pour des oi^ganes de 
regulation thermique de 1'appareillage de synthase. II sera apprecie que 
cette fagon de proceder permet d'intervenir sur les organes de regulation 
15 thermique de 1'appareillage de synthase avant meme que la temperature ne 
change. Des mesures de temperatures interviennent par ailleurs 
avantageusement comme grandeurs d'entree suppiementaires dans 
ralgorithme esclave. 

L'algorithme maitre comprend avantageusement la structure 
20 suivante : 

- un organe de prediction, base sur un module direct du precede permettant 
de fournir une prediction des proprietes du produit synthetise en fonction de 
mesures et/ou de predictions des concentrations des constituants ; 

- un organe d'adaptation comparant les predictions de proprietes calcuiees 
25 par 1'organe de prediction a des valeurs effectivement mesurees sur le 

produit synthetise et deduisant de cette comparaison des parametres 
d'adaptation, lesdits paramfctres d'adaptation intervenant comme grandeurs 
d'entree suppiementaires dans ledit organe de prediction de l'algorithme 
maitre ; et 

30 - un organe de controle, base sur un module inverse du procede, pour 
calculer, en fonction des consignes et des predictions de proprietes du 
produit k synthetiser, des consignes de concentrations pour Falgorithme 
esclave, lesdits parametres d'adaptation intervenant egalement comme 
grandeurs d'entree suppiementaires dans ledit organe de contrSle. 

35 L'algorithme esclave comprend avantageusement la structure 

suivante : 
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- un organe de prediction, base sur un modele direct du precede permettant 
de foumir une prediction des concentrations d'un ou plusieurs des 
constituants sur base d'un bilan de matiere dans le reacteur ; 

- un organe d'adaptation, comparant les predictions de concentrations 
calculees par le modele direct a des mesures de concentration et deduisant 
de cette comparaison des parametres d'adaptation, Iesdits parametres 
d'adaptation intervenant comme grandeurs d'entree supplemental dans 
ledit oigane de prediction de 1'algorithme esclave ; et 

- un oigane de controle, base sur un modele inverse du precede, pour 
calculer, en fonction de la consigne d'allure de marche, des consignee de 
concentration calculees par l'organe de controle de ralgorithme maitre et 
des predictions de concentration calculees par l'organe de prediction de 
ralgorithme esclave, les consignes pour les debits entrant dans le 
reacteur, Iesdits parametres d'adaptation intervenant comme grandeurs 

15 d'entree supplementaires dans ledit organe de controle de I'algorithme 

esclave. . 

La dynamique du precede est avantageusement decrite et calculee au 
moyen de fonctions du type y = LAG (u , » 

cette fonction etant la solution de l'equation differentielle 



10 



20 



U = T • + y 

dt * 



ou les arguments u et t varient dans le temps. L'utilisation de cette 
fonction conformement aux theoremes 1 et 2 exposes ci-dessous peimet de 
resoudre de maniere sequentielle les bilans de matiere utilises par 
ralgorithme esclave et de decrire la cinetique du precede par des equations 
caracteristiques simples dans ralgorithme maitre. La fonction LAG permet 
25 encore de neduire considerablement le volume des calculs necessaires et 
rend par consequent inutile l'utilisation d'ordinateurs rapides et puissants. 
De plus, cette fonction permet d'etablir de facon particuUerement simple les 
modeles directs et inverses du precede ou de certaines de ses parties. 
Les principals qualitis de la regulation proposee peuvent etre 
30 resumees comme suit : 

- anticipation : la regulation commence a corriger les perturbations 
mesurees avant meme que leur effet ne se soil manifeste sur les mesures 
■ des proprietes (utilisation de predictions de concentrations, de predictions 
de proprietes et de predictions de temperatures dans les algorithmes); 



- precision meme en presence de perturbations : le module direct et le 
module inverse sont recalibres en permanence en utilisant les mesures des 
proprietes (adaptation); 

- validity etendue : 1'algorithroe garde sa validity durant les transitions 
d'allure de marche et de qualite, ainsi que lors des demarrages et des 
arrets (raise en equation de la dynamique du proc&te, utilisation de 
predictions pour les grandeurs dont les mesures impliquent des temps 
moits importants) ; 

- simplicity : le d^veloppement et la mise en oeuvre sont facility grace a 
une m&hode originate de mise en Equation de la dynamique du precede 
(fonction LAG). 

Le proc&te de synthese a nSguler est ainsi modyiise sous une fonne 
qui est gendralement qualifiee de "modfcle de connaissance" (Tim principle 
model"), c'est-i-dire que son module est dlabord a partir d'&juations 
refletant le dSroulement physico-chimique dftaillg du precede. Une telle 
approche pennet d'obtenir, au moyen d'un jeu d'equations relativement 
simple sur le plan math&natique, des n&ultats sup^rieurs a ceux qui auraient 
6t6 obtenus au moyen d'un modfcle empirique de type boite noire, en 
foumissant notamment des paramfetres relies h des grandeurs i^elles et une 
meilleure validity k l'exterieur de Tespace d'identification (extrapolation). 
La plupart des modeles empiriques utilisent des equations complexes, 
souvent d'ordre eleve si Ton desire obtenir une simulation correcte de la 
dynamique du precede, et dont les paramfetres (notamment les constantes de 
temps) doivent 6tre identifies pour un point de fonctibnnement precis; le 
modele n'est valable que dans le voisinage imm&iiat de ce point de 
fonctionnement. Une telle approche peut difficilement etre generaiisee h un 
grand nombre de points de fonctionnement dans le cas d'un precede de 
synthfese chimique reel faisant intervenir de nombreuses grandeurs. 

Au contraire, selon la methode de regulation de r invention, on 
utilise un ensemble d'equations simples, purement statiques ; la dynamique 
du proc&te est simulee par des fonctions simples (cf . la fonction LAG ci- 
dessus). Avantageusement, les temps de sejour (constantes de temps des 
Equations) peuvent etre recalcutes aussi souvent qu'on le souhaite, ce qui ne 
pose aucun probfeme etant donne la simplicity des equations. On obtient en 
fin de compte un ensemble d'equations extremement simples et faciles a 
resoudre en temps reel, meme k frequence yievee. 
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La methode de regulation propose* est avantageusement applicable a 
des process de synthese, notamment en continu, de polymeres 
(polymerisation), et notamment a la polymerisation en continu d'olefines 
teUes que par exemple r^thylene ou Ie propylene, aussi bien en phase 
liquide qu'en phase gazeuse. 

La presente invention concerne egalement un precede de synthese 
d'un ou plusieurs produits chimiques, reguie au moyen d'une methode de 
regulation conforme a l'invention. En particulier, de tits bons nfcultats ont 
ete obtenus pour la regulation d'un proceed de synthese en continu de 
polyethylene par polymerisation d'ethylene dans au moins un reacteur les 
reactifs comprenant Tethylene, 1'hydregene et/ou un comonomere optionnel 
la reaction de polymerisation ayant lieu dans un solvant en presence d'un 
catalyseur et une partie du contenu du reacteur etant proleve en permanence 
ou par intermittence. Ce precede pent aussi bien se derouler en phase 
liquide qu'en phase gazeuse ; de preference, il se deroule en phase Uquide 
(dans un solvant). 

La methode s'applique de maniere analogue a la synthese du 
polypropylene (le monomere de depart principal etant dans ce cas le 
propylene au lieu de F&hylene), du propane pouvant element etre present 
s. le procede se deroule en phase gazeuse. Pour le polypropylene, 1'indice 
de Hmdite en fondu est souvent designe par MFI au lieu de MI. 

L'invention concerne egalement un dispositif de regulation destine a 
mettre en oeuvre la methode de regulation de l'invention, ainsi qu'une 
installation de synthese d'un ou plusieurs produits chimiques comprenant un 
25 tel dispositif de regulation. 

De maniere plus precise, l'invention concerne egalement un 
dispositif de regulation d'un procede de synthese d'un produit chimique dans 
un appareillage de synthese comprenant au moins un reacteur, ledit 
dispositif comprenant : 
30 - au moins une unite de calcul; 

- desmoyens pour entrer des consignes de propria du produit a synthe- 
tiser dans l'unite de calcul; 
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- des moyens pour entrer des consignes d'allure de marche du produit a 
synthetiser dans 1' unite de calcul; 

- des organes de mesure de debit des flux entrant dans le reacteur 

- des organes de mesure de la composition des flux sortant du reacteur; 
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- des organes de commande de debits (actionneurs) pour regler les debits 
des flux entrant dans le reacteur; 

- des moyens de communication entre ladite unite de calcul, lesdits organes 
de mesures de debit et lesdits oiganes de regulation; 

dans lequel : 

- la masse d'au moins un constituant est calcuiee par la fonction LAG 
comme expose ci-dessus; 

- 1'unite de calcul est apte a calculer a l'aide d'un algorithme maitre, en 
fonction des consignes de propriet^s, des consignes de concentration de 
reactifs dans le reacteur; 

- 1'unite de calcul est apte k calculer k Taide d'un algorithme esclave, en 
fonction des consignes de production et des consignes de concentration, 
de consignes de debit pour les flux entrant dans le reacteur, ces 
consignes de debit etant transmises comme consignes cT entree aux 
organes de regulation de debits; 

- les mesures r&lis&s par les organes de mesure de debits inteiviennent 
comme grandeurs d'entree supplementaires dans ledit algorithme esclave, 
pour permettre a ce dernier de calculer, en fonction de ces mesures de 
debits, des predictions de concentrations, ces predictions de 
concentrations intervenant dans V algorithme maitre pour calculer des 
predictions de propriet6s utilisees comme grandeurs d'entree 
supplementaires dans le calcul des consignes de concentrations. 

L' invention concerne egalement un dispositif de regulation tel que 
d&rit ci-dessus, dans lequel : 

- I'appareillage de synthese comprend en outre : 

des organes de regulation thermique apte a contr61er la temperature 

dans le reacteur, et 

des capteurs de temperatures; 

- F algorithme maitre est apte a calculer, en fonction des consignes de 
propri&es, des consignes de temperatures pour le reacteur; 

- Talgorithme esclave est apte a : 

calculer un bilan thermique pour le reacteur, 
- resoudre ce ou ces bilans thermiques de fa?on a determiner la 

chaleur qu'il faut ajouter a, respectivement retrancher de la synthese 

pour respecter les consignes de temperature, et 

deduire de ce ou ces bilans thermiques des consignes d'entree pour 
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les oiganes de regulation thennique du reacteur, et 
- recevoir comme grandeurs d'entree supplementaires les mesures 
effectuees par les capteurs de temperature. 

L'invention conceme egalement un dispositif tel que decrit ci-dessus, 
dans lequel des mesures effectuees par les capteurs de temperature 
interviennent comme grandeurs d'entree supplementaires dans ralgorithme 
maitre. 

L'invention conceme egalement un dispositif tel que decrit ci-dessus, 
comprenant : 

- au moins un analyseur apte a foumir des mesures des proprietes inter- 
venant dans 1'algorithme maitre; et 

- des moyens pour entrer ces mesures de propriety dans 1'unite de calcul; 
ladite unite de calcul comprenant : 

- un premier organe de prediction, base sur un premier modele direct du 
precede pennettant la prediction des proprides du produit synthetise, en 
fonction des predictions de concentrations calculees par ralgorithme ' 
esclave; 

- un premier organe d'adaptation, comparant les predictions de proprietes 
calculees par le premier organe de prediction aux valeurs effectivement 

20 mesurees sur le produit synthetise et detiuisant de cette comparaison des 

parametres d'adaptation, intervenant comme grandeurs d'entree 
supplementaires dans ledit premier organe de prediction ; et 

- un premier organe de conrrole, base sur un premier modele inverse, pour 
calculer, en fonction des consignes et des predictions de proprietes, des 

25 consignes de concentrations pour ralgorithme esclave, lesdits parametres 

d'adaptation intervenant egalement comme grandeurs d'entree 
supplementaires dans ledit premier organe de conrrole. 

L'invention conceme egalement un dispositif tel que decrit ci-dessus, 
comprenant en outre : 

30 - au moins un analyseur apte a foumir des mesures de concentration des 
reactifs; et 

- des moyens pour entrer ces mesures de concentration dans I'unke de 
calcul; 

ladite unite 1 de calcul comprenant : 
35 - un deuxieme organe de prediction, base sur un deuxieme modele direct 
pennettant la prediction des concentrations, en fonction du bilan de 
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matifere dans le reacteur; 

- un deuxieme organe d'adaptation, comparant les predictions de 
concentrations calculees par le deuxieme organe de prediction aux 
mesures de concentration et d<Sduisant de cette comparison d'autres 

5 parametres d'adaptation, intervenant comme grandeurs d'entree 

supplementaires dans ledit deuxieme organe de prediction; et 

- un deuxifeme organe de controle, bas£ sur un deuxieme modfele inverse, 
pour calculer, en fonction des consignes de production, des consignes de 
concentration calculees par 1'algorithme maitre et des predictions de 

10 concentration du deuxieme organe de prediction, des consignes pour les 

debits entrant dans le reacteur, lesdits autres parametres d'adaptation 
intervenant comme grandeurs d'entree supplementaires dans ledit 
deuxieme organe de controle. 

Un mode de realisation concret de 1' invention est illustre sur base 
15 d'un procedd de synthese en continu de polyethylene (PE), en se reterant 
aux Figures 1 k 10. Celles-ci montrent : 

Figure 1 : un schema d'un circuit de fabrication de polyethylene; 
Figure 2 : un schema simplify de la structure d'une regulation 
avancde selon 1'invention; 
20 Figure 3 : un schema de principe de la regulation avancde appliquee 

au circuit de fabrication de la Figure 1; 
Figure 4 un schema de principe d'un algorithme de controle 

adaptatif, tel qu'utilis6 dans le systeme de regulation 
avanc£e selon la Figure 2; 
25 Figure 5 un schema de la structure de 1'algorithme maitre dans le 

systeme de regulation avancde selon la Figure 2; 
Figure 6 : un schema de la structure de 1'algorithme esclave dans le 

systfeme de regulation avanc& selon la Figure 2; 
Figure 7 : le schema general d'une methode de regulation conforme 
30 k 1'invention; 

Fig. 8-10 : le schema de variantes particulieres de la methode de 
1'invention. 

Sur la figure 7, on distingue tout d'abord le proc&ie de synthase 
proprement dit (Pr), qui peut etre contrdte en imposant au moins une 
35 consigne d'une grandeur de commande (CGC) (par exemple un ou plusieurs 
debits de constituants entrant dans le reacteur, une temperature, etc.) k un 
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actionneur adequat (vanhe, dispositif de chauffage ou de refroidissement, 
etc.). La regulation s'effectue au moyen d'un organe de controle (OG), 
base" sur le modele mathematique inverse du precede, et dont le role 
principal est de comparer la ou les consignes des grandeurs reglees (CGR) 
(par exemple une ou plusieurs proprietes du produit a synthetiser et/ou une 
ou plusieurs grandeurs liees au ddroulement du precede) a la ou aux 
predictions de ces grandeurs (PGR). La ou les predictions des grandeurs 
reglees (PGR) sont calculees par un organe de prediction (OP), base sur le 
modele mathematique direct du precede, sur la base de mesures des 
grandeurs de commande (MGQ. On notera qu'aucune mesure de 
propriete(s) du produit synthetise n'intervient dans la regulation. 

La figure 8 represente une variante de la methode de la figure 7, 
dans laqueUe le modele mathematique du precede est periodiquement adapte 
par un organe d'adaptation (OA), sur la base de Pecan (avantageuseraent 
15 fdtre, ou traite numeriquement) entre les predictions (PGR) et les mesures 
(MGR) des grandeurs reglees. Une resynchromsation (recalage dans le 
temps) de ces mesures et de ces predictions est souvent necessaire, par 
exemple lorsqu'il s'agit de mesures de proprietes du produit synthetise dont 
la duree d'obtention est longue. L'organe d'adaptation (OA) transmet le 
20 resultat de ses calculs, c'est-a-dire ses instructions d'adaptation, a l'organe 
de prediction (pour adapter le modele direct du precede) et a l'organe de 
controle (pour adapter le modele inverse du precede). On notera que la ou 
les mesures de proprieteXs) du produit synthetise ne sont pris(es) en compte 
que dans le processus d'adaptation, qui se deroule generalement a une 
25 frequence beaucoup plus foible que le processus normal de controle. La 
lenteur eventueUe de ces mesures n'a done aucun effet direct sur la qualite 
de la regulation. 

La figure 9 represente une autre variante de l'invention, dans 
laqueUe une ou plusieurs mesures de grandeurs reglees (Mq R ) sont prises 

30 en compte par l'organe de controle (OQ, et une ou plusieurs mesures de 
grandeurs reglees (M GR *) (eventuellement differentes) sont prises en 
compte par l'organe de prediction (OP). De meme, une ou plusieurs 
mesures de grandeurs de commande (M QC ') peuvent en outre etre prises en 
compte par l'organe de controle (OC). 

35 11 v » soi qu'on poumit creer une autre variante de la methode de 

l'invention en combinant les variantes des figures 8 et 9. c'est-a-dire a la 
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fois en utilisant un organe d'adaptation et en prcnant en compte une ou 
plusieurs mesures de grandeurs de commande dans 1'organe de contrdle, 
et/ou une ou plusieurs grandeurs de commande dans 1'organe de prediction 
et/ou dans 1'organe de contrive. 
5 Dans la figure 10, le module math&natique du procede n'est pas 

adapte au sens propre du terme, mais on utilise T&art (avantageusement 
filtre) entrc les mesures et les predictions des grandeurs reglees pour 
corriger les consignes des grandeurs rfgl6es (C GR ). En l'occurrence, cette 
correction a &6 ici representee comme une simple difference : on soustiait k 
10 chacune des consignes de grandeurs regl&s un terme coirectif calculd par 
1'organe d'adaptation OA (qui n'est en fait ici qu'un organe de correction), 
ce qui foumit des consignes corrigees C GR \ transmises a J'oigane de 
controle OC. H va de soi que dans certains cas, la correction peut 
impliquer des operations plus complexes qu'une soustraction, par exemple 
15 une division (dans ce cas-ci, on peut cependant se ramener a une 

soustraction en considerant les logarithmes des grandeurs considgrees). 
Cette methode est communement appelfe Internal Model Control (IMC). 

Si Ton se rgfere k la Figure 1, qui represente schematiquement un 
circuit de synthase en continu de polyethylene (PE), la polymerisation de 
20 rethylene se fait dans un reacteur boucle 10, en suspension dans un solvant 
ad^quat, tel que par exemple l'hexane. Le processus est continu, c'est-a- 
dire que les reactifs sont injectes en continu, et qu'une partie du contenu du 
reacteur 10 est continuellement soutirfe. Une pompe de circulation (non 
representee) assure Phomogeneite du contenu du r6acteur 10. 
25 Les reactifs introduits dans le reacteur sont Fethylene "Et w , 

Thydrogene tt Hy" et le butfene "Bt" (cf. reference 11). Un catalyseur est 
egalement inject^ en continu. H importe de bien maitriser les concentrations 
des reactifs dans le reacteur, car les proprietes de la resine PE sont 
principalement d6terminees par les rapports de concentrations Hy/Et et 
30 Bt/Et. 

La temperature de polymerisation dans le reacteur est un paramfetre 
additionnel qui influence les proprietes de la resine PE. Vu que la reaction 
de polymerisation est fortement exothermique, la temperature du reacteur 
doit etre rgguiee en ayant recours k un circuit de refroidissement 12. 
35 ^ reacteur 10 en fonctionnement contient done du solvant, du poly- 

mere et des reactifs n'ayant pas encore reagi et du catalyseur. Son contenu 
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est soutire en permanence a travers la conduite de soutirage 14. Ce contenu 
soutire entre dans un stripper STP 16, qui separe le polymere PE et les 
fluides (solvant et reactifs). Ces fluides sont vaporises par injection de 
vapeur d'eau et evacues dans un condenseur CD 18. Dans ce dernier le 
solvant est de nouveau condense avant d'etre recycle. Les reactifs qui sont 
plus legers sont separes du solvant et sont aussi recycled Un chromato- 
graphe en phase gazeuse GC (20) place a la sortie du condenseur 18 permet 
de determiner les concentrations Hy/Et, Bt/Et des reactifs. 

Le polymere evacue du stripper 16 est concentre dans une 
centrifugeuse CFG 22, puis seche dans un secheur a lit fluide SHLF 24 
avant d'etre envoye vers le finishing pour y etre granule. Des echantillons 
sont preleves a la sortie du secheur 24 pour y mesurer les proprietes de la 
nJsine : cristallinite (mesuree par la masse volumique specifique "MVS") et 
proprietes rheologiques (indice de fluidite en fondu ("melt index (MI)" ou 
"melt flow index (MFD") et viscosite en fondu > 2 \ mesuree dans un 
ecoulement capillaire sous un cisaillement de 100 s" 1 ). 

Ia dyiiamique de ce processus de synthese de PE est lente et 
complexe : 

- Le reacteur boucle 10 se comporte comme un reacteur parfaitement 
melange. Des lore, toute modification du detit d'alimentation d'un des 
reactifs ne se traduira que progressivement sur la concentration de ce 
reactif dans le reacteur. En effet, le nouveau d6bit doit se melanger au 
volume tout entier du reaaeur 10 pour Tamener a la nouveUe concen- 
tration d'6quilibrc. 

?. La mesure des concentrations des reactifs se fait par un chromatographe 
en phase gazeuse 20; U s'agit d'un appareil discontinu, qui procede par 
etapes successives : prelevement d'un echantillon de gar, analyse, puis 
elaboration des resultats. D y a done un temps mort (de 5 a 15 niin.) 
entre les changements de concentration et leur mesure. 

- Les proprietes du polymere fabrique a chaque instant dependent piincipa- 
lement des concentrations des reactifs. Toute modification de ces 
concentrations affecte done instantanement les proprietes du polymere 
fabrique\ Par contre, les proprietes moyennes dans le reacteur ne se 
modifient que progressivement, car le polymere fraichement produit doit 
se melanger au polymere deja present dans le reacteur 10 (temps de 
sejour : ± 2 h). 
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- Lorsque le polymfere est soutirg du reacteur 10, il subit a nouveau une 
serie de melanges dans les differents appareils (STP, CFG et SELF) 
•destines a le s&her (temps de s6jouf : ± 2 h). Des echantillons de 
polymfere sont alors pr^leves et analyses par le laboratoire de l'usine. 
5 . Les r&ultats de ces analyses ne seront done communiques k la fabrication 
qu'apres un nouveau temps moit, qui peut etre important (± 2 h). 
Mod^isation du precede k l'aide d'une fonction LAG 

Selon la methode de l'invention, la modelisation dynamique d*un 
proc&te de synthase en continu est realis6e en ayant recours k des hypo- 
10 theses de melanges parfaits et de d61ais purs. Les melanges parfaits sont 
mis en Equation au moyen d'une fonction bien connue des ing&iieure, la 
fonction "LAG", ou filtre passe-bas (du ler ordre) ; cette fonction est 
lindaire et aisement programmable. 
EUe est d6finie comme suit : y = LAG (u , r) 
15 (dit "LAG de u durant r") comme etant la solution de r equation 
difKrentielle 

u = T • + y 
dt 

dont les arguments u et r varient dans le temps. 

Cette Equation peut se resoudre num&iquement (meme en temps 
r6el) par une equation alggbrique du ler ordre, ayant pour arguments les 
20 variables suivantes : 

- la periode d'&hantillonnage M T" (ou temps ecoule depuis la derniere 
iteration) 

- le temps de sejour (ou "constante" de temps) "r" a Tinstant "t" 

- la grandeur d'&at "y" a Tinstant precedent "t-T" 
25 - la grandeur de commande "u" a l'instant present "t" 

(u et r represented en fait les valeurs mesurees ou calculees k r instant "t" 
sur les grandeurs V et "t w , que Ton suppose avoir 6te constantes dans tout 
TintervaUe precedent "T). 

De preference, T est petit par rapport a r (par exemple au moins 10 
30 fois plus petit), de fagon a augmenter la precision du calcul. 

La solution de F equation susmentionnee peut se calculer par exemple 
par la formule suivante : 

y(t) = y(t-T) . e - T/r W + u(t) . ( 1 - e ' T,T ® ) 
ou encore, plus simplement (approximativement) : 
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y(t-T) + u('t)-— I- 
y(t') = T (t) 



La modelisation du precede a l'aide de la fonction LAG se base sur 
Ies theoremes suivants : 

Soit un reacteur parfaitement melange" (CSTR), de volume V R . 
Divers constituants (reactifs ou inertes) ralimentent, paraii lesquels le 
reactif "x" (deTrit entrant Fx m ) ayant la propriety "Px^" a I'entree. On 
mesure egalement le d£bit de sortie "Four" (soutirage). 
Thforftme 1 : Application de la methode LAG au calcul d'un bilan 
massique : 

A chaque instant, la masse "Mx R " d'un constituant "x" dans un 
reacteur parfaitement melange (CSTR) est egale au LAG du produit du deM 
massique entrant "Fx^" par un temps V x ", durant ce meme temps "r x " : 
Mx R = LAG(Fx IN . r x ,r x ) ' (en kg) 

Le temps V x " est le "temps de sejour de x"; il vaut la masse du constituant 
Mx R , divisee par la somme des debits massiques "sonants" (quantite cons- 
15 ommee par reaction "FRx", debit sonant du reacteur "Fxqjjj.", etc.). 

r x « Mx R / (Fx 0UT + FRx + ...) ' * (en h) 

Ce theoreme donne ainsi une methode exacte pour le calcul dynamique 
(meme en temps reel) des concentrations dans un reacteur parfaitement 
melange. En effet, la concentration "Cx R " du constituant "x" exprimee en 
20 kg/m 3 vaut, si **V R " est le volume du reacteur : 

Cx R = Mx R /V R (enkg/m 3 ) 
Connaissant le d&nt volumique total "FV OUT " sonant du reacteur, 
I on definit le temps de sejour du reacteur "t r " : 
r R = V R / FV OUT 

25 Des lors, le debit massique "Fxquj" du constituant "x" en sortie du 
reacteur (soutirage) vaut : 

Fx OUT = Mx R /7 R 

On remaiquera que si "x" est un inerte (ne subissant pas de reaction 
et ne sonant du reacteur que par soutirage), on a : 
30 r x = r R 

Par ailleurs, dans le cas frequent oil la vitesse de reaction de "x" est 
proponionnelle a sa concentration Cx R) avec un facteur de propoitionnalite 
R x> on a : FRx = R x . Mx R et dbs lors 
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r x = 1 / (R X + l/r R ) 
Thforfeme 2 : Application de la methode LAG au calcul d'une propriete 
d'un mdiange : 

Soit une propria M Px w d'un constituant, repondant a la loi de 
5 m61ange lindaire : 

Px 1+2 = wj . Pxj 4- w 2 . Px 2 
oii W| et w 2 sont les fractions massiques de propria Pxj et Px 2v (avec Wj 
+ w 2 = 1). 

A chaque instant la propriety Px R dans un reacteur parfaitement 
10 melange (CSTR) est 6gale au LAG de la propriety a l'entrfe Px^, durant 
un temps de s^jour qui vaut le rapport de la masse Mx R du constituant dans 
le r&cteur divis6 par le debit massique entrant (et/ou apparaissant par 
reaction) Fx^ : 

Px R = LAG ( Pxjn , Mx R / Fx m ) 
15 II est ainsi possible de tenir compte de la dynamique du precede, et 

de recalculer en permanence ses constantes de temps. 

Comme expos6 prec^demment, les "prop^i&& ,, dont il est question 

ici peuvent dans certains cas avoir subi une transformation math^matique les 

rendant lineaires (par exemple, le logarithme de 1'indice de fluidity en fondu 
20 d'un polymere peut etre considere comme repondant k une loi de melange 

lin£aire). 

Principe de la regulation : 

Le module du proc&te de synthase etant &abli, il faut un algorithme 
qui calcule les parametres necessaires k la regulation de ce proc&te. 

25 La Figure 2 montre un schema general simplify du type de 

regulation avancfe (Advanced Process Control ou "APC") adopts pour le 
proc&ie de polymerisation decrit plus haut. On voit que ce syst&me de 
regulation comprend une cascade de deux algorithmes, cette cascade pilotant 
notamment les regulateurs PID des debits d'alimentation en reactifs. 

30 Les deux algorithmes en cascade, appeles algorithme maitre et algo- 

rithme esclave, sont tous deux des algorithmes dynamiques adaptatifs bases 
sur des modeles issus de la connaissance du processus (par opposition aux 
modeles empiriques), bases notamment sur les bilans de mat&re et la 
cinetique du precede regule. lis utilisent la fonction LAG predefinie. 

35 Sur la Figure 3, qui montre le principe du systfeme de regulation 

dans le contexte du precede de polymerisation d6crit plus haut, on voit 
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que : 

- I'algorithme maitre est base sur les equations caracteristiques des cata- 
lyseurs, c'est-a-dire des equations foumissant les proprietes du PE en 
fonction de la temperature de polymerisation et des concentrations des 

5 reactifs dans le nSacteur; il fournit a I'algorithme esclave les consignes 

des concentrations en reactifs pour satisfaire les consignes pour les 
.propriety du PE; 

- Talgorithme esclave est base sur un bilan matifere et la cinetique chimique 
des reactions; il fournit aux regulateurs des debits d'alimentation les 

10 consignes de debits de reactifs n&essaires pour satisfaire les consignes de 

concentrations impos&s par Falgorithme maitre et la consigne d'aliure de 
marche du precede. De preference, il calcule egalement un terme 
anticipatif ( tt feed-forwanT) pour la consigne de temperature, ameliorant 
la regulation de temperature, lors des changements d'allure de marche. 
15 Ce type de regulation n'est d'une precision parfaite que si le modele 

est parfait et prend en compte toute les perturbations possibles. Cela n'est 
en general pas le cas. C'est pouiquoi en general (voir Figure 4) on adapte 
en continu le modele direct (et le modfele inverse) en comparant les 
predictions aux mesures des proprietes. Cette "adaptation" du modfele 
20 permet de maintenir sa precision en presence de perturbations non 
modelisees, et ainsi d'obtenir une regulation plus precise en toutes 
circonstances. 
L'algorithme esclave 

La Figure 6 illustre le principe de ralgorithme esclave : 
25 1. un organe de prediction utilisant un modele direct du precede prddit, a 
partir des mesures des debits alimentant le reacteur, les concentrations 
en reactifs et en polymere; 

2. un organe d'adaptation compare les concentrations en ethylene (Et), 
hydrogene (Hy) et butene (Bt) mesuree par un analyseur (chroma- 

30 tographe k gaz) aux valeurs pnSdites par le modele direct, de fxpn k 

determiner trois parametres d'adaptation : 

- Tactivite specifique du catalyseur pour r ethylene "AEt", en kg/h de 
polyethylene par kg de catalyseur et par kg/m 3 d'ethylene 
Terreur de gain sur la mesure de debit d'hydrogene "KfHy" 
35 - la purete de Talimentation en butfene "KfBt" 

3. Torgane de controle calcule; a panir des consignes de concentration 
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calculees par l'algorithme maitre et de la consigne d'allure de marche, 
les consignes pour les debits d'alimentation des reactifs ; ces consignes 
de ddbit sont composees d'un feed-forward base sur le modfcle inverse, 
et d'un feed-back proportionnel a I'ecart entre le module direct et les 
5 . consignes de concentration. 

Pour comprendre les calculs effectu& par l'algorithme esclave, dans 
le cas d'un proc&te de synthase de polyethylene, il feut d'aboid savoir qu'il 
est generalement admis que la vitesse de polymerisation "VitPj^" est 
proportionnelle : 

10 - I la concentration en Et non poIym6ris<5 cEt R (en kg/m 3 ), 

- k la concentration en catalyseur actif dans le reacteur cCA R (en kg/m 3 ), 
et 

- a un facteur de proportionnalitd, fonction (difficilement quantifiable) de 
la temperature, des concentrations en Hy, en Bt et en co-catalyseur, de la 

15 pr6sence de poisons, etc. Ce tenne est appeld "activity catalytique" pour 

T&hylfene Aj^. Sauf dysfonctionnement majeur (poison, etc.), il varie 
raisonnablement peu en cours de campagne. 
VitP a = A B . cCA R . cEt R (kg/m 3 .h) 

La quantite d'&hylfene TpR" qui polymerise dans le nSacteur a 
20 chaque unitg de temps ( M d6bit de polymerisation") vaut done, si V R est le 
volume du reacteur : 

FpEt = VitP^. V R (kg/h) 

= Agj . cCA R . MEt R 
ou MEt R est la masse d'&hylfene en solution dans le reacteur (en kg) 
25 On sait par ailleurs que la vitesse d'incoiporation de l'Hy est environ 

100 fois plus lente que 1'Et, et 10 fdis plus lente pour le Bt. On en deduit : 
FpHy = A m . cCA R . MHy R / 100 
FpBt = A a . cCA R . MBt R / 10 

ou MHy R est la masse d'hydrogfcne en solution dans le reacteur (en kg) et 
30 MBt R est la masse de butfene en solution dans le reacteur (en kg) 

Le modele esclave utilise maintenant les mesures suivantes : 
FEtjjj = debit d'alimentation en ethylene (monomfere) (kg/h) 
FSv IN = debit d'alimentation en solvant (hexane) (kg/h) 
FCA IN = d6bit d'alimentation en catalyseur (kg/h) 
35 FHyjj^ debit d'alimentation en hydrogfene (agent de transfert) (kg/h) 
FBt^ = d6bit d'alimentation en butfene (co-monomfere) (kg/h) 
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II utilise egalement les parametres d'adaptation suivants : 
Agj = "activity catalytique" pour fethylene 
KfHy = erreur de gain sur la mesure de ['alimentation en hydrogene 
KfBt = puret£ de l'alimentation en butene. 
5 . Les calculs suivants sont effectues de facon sequentielle, dans Foidre 
suivant, avec une frequence elevfe (l'intervalle de temps separant chaque 
iteration etant faible par rapport au plus petit des temps de sejour r x ). 

Le volume du reacteur etant constant, le debit volumique soitant 
equivaut au debit volumique entrant (fluides.incompressibles). Le debit 
10 volumique sortant "FV 0UT " peut de* lore etre calcule comme la somme 

des debits massiques entrants divises par la masse volumique qu'ils ont dans 
le reacteur : 

FV OUT = FS VlN /650 + FEt m /950 + FBt^/600 ( m 3/ h) 
(ou les masses volumiques sont les suivantes : 650 kg/m 3 pour le solvant, 

15 950 kg/m 3 pour le polyethylene, 600 kg/m 3 pour le butene). On fait ici 
l'hypothese que tout T&hylene se transforme instantanement en poly- 
ethylene, et on neglige le debit d'hydrogene et de catalyseur (quelques kg). 

Le solvant est chimiquement inerte et ne sort du reacteur que par 
soutirage du reacteur. On calcule sa masse "MSv R " dans le reacteur en 

20 utilisant le theoribme 1 : 

T R = V R ' ^OUT 00 (temps de sejour dans le reacteur) 

MSv R = LAG (FS VlN • r R>TR ) (kg) 

Le catalyseur se desactive avec une constante de temps "kd" ; la 
masse "MCA R " de catalyseur actif dans le reacteur est calculee comme 
25 suit : 

r CA " 1 1 Wtr + kd) (h) 
MCA R = LAG (FCAjn . r CA , r CA ) (kg) 
et done la concentration "cCA R " de catalyseur actif dans le reacteur vaut ■ 
cCA R = MCA R /V R (kg/m 3 ) 
30 Methylene "sort" du reacteur dans le debit de soutirage ainsi que 

dans la reaction de polymerisation. Sa masse dans le reacteur "MEt R " est 
calculee comme suit : 

t r = 1 / (l/r R + A a . cGA R ) (h) 
MEt R = LAG (FEtjN . , r a ) (kg) 
35 D'une racon similaire, on calcule la masse "brute" (non calibree) 

d'hydrogfene dans le reacteur "MHy RAW " : 
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r Hy = 1 / (l/r R + A B . cCA R / 100) (h) 

MHyRAW = 1LAG (F^IN * T Hy - T Hy ) (*8> 
La masse "MHy R ", corrigee pour prendre en compte Ferreur de 

gain sur la mesure de Talimentation en hydrogene, vaut : 

5 MHy R = KfHy . MHy RAW (kg) 

Le rapport Hy:Et dans le reacteur vaut done : 

HyEt R = MHy R / MEt R 

D'une fa?on similaire, on calcule la masse de butene "bnite" 

,, MBt RAW w : 

10 r Bt = 1 / (1/t r + A a . cCA R / 10) (h) 
MBt RAW = LAG (FBt m . r Bt , r Bt ) (kg) 

La masse tt MBt R ", corrigee pour prendre en compte la purete de 
I'alimentation en bufene, vaut : 

MBt R = Kffit . MBt RAW (kg) 
15 Le rapport Bt:Et dans le reacteur vaut done : 

BtEt R = MBt R / MEt R 

On a vu que le dibit de polymerisation TpR" (allure de marehe 
instantanfe) vaut : 

FpEt = A a . cCA R . MEt R (kg/h) 
20 Le polymere etant inerte et ne subissant pas de reaction, sa masse 

dans le reacteur "MPE R " vaut : 

MPE R = LAG ( FpEt . t r , r R ) (kg) 
Le debit de polymere sortant du reacteur "FPEqut" vaut des lors : 
FPEquj = MPE r / r R (kg/h) 

25 Adaptation du modfele esclave : 

Le bloc d'adaptation utilise un analyseur (par exemple un chroma- 
tographe & gaz "GC") pour obtenir les mesures des concentrations dans le 
reacteur en ethylene "cBLqq , en hydrogene "cHy^" et en butene 
McBt GC" (exprim&s par exemple en kg/m 3 ). Ces mesures sont compares 

30 aux valeurs predites par le module direct, afin de determiner les trois 
paramfetres d'adaptation suivants : 

- Tactivite sp^cifique du catalyseur pour l'ethyfene "AEt n , en kg/h de 
polyethylene par kg de catalyseur et par kg/m 3 d'ethylfene 

- Terreur de gain sur la mesure de debit d'hydrogene "KfHy" 
35 - la purete de I'alimentation en butfene "KfBt" 

Le chromatographe a gaz fournit des mesures echantillonnees avec 
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un retard d'environ 6 minutes. 

Le calcul de I'activite sp&ifique pour l'&hylene "AEt" d&oule des 
egalit& suivantes : 

cEt GC =cEt R (kg/n,3) 
= MEt R /V R 

= LAG <FEt m • 7^ , r a ) / V R 

= LAG (FEt^ / (l/r R + A a . cCA R ) , r a ) / V R 
- A m " = ( FEt^ / LEAD ( cEt^ . V R , r a ) - l/r R ) / cCA R 

La mesure cEt^ etant echantillonnee et bruitee, on ne souhaite pas 
la trouver dans un LEAD; des lors, on utilise plutdt : 
V LAG (A E - , r B ) 

■ . = LAG (FEtjj^/cCA R , r a ) / (cEt Gc .V R ) - 1 / LAG (r R .cCA R , r E ) 

On tient compte du retard de 6 minutes sur la mesure de cEt GC en 
introduisant deux LAGs en serie de 3 minutes chacun sur les valeure du 
15 modele, et on obtient la formule finale : ■ 

A e( = LAG ( FEt m / cCA R , r a , 3/60 , 3/60 ) / (cEt^ . V R ) 
- 1 / LAG ( t r . cCA R ,T m , 3/60 , 3/60) 
• Le calcul du gain du dfbit d'hydrogene "Kffly" decoule des egalites 
suivantes : 

20 cHy GC =cHy R (kg/ m 3) 
= MHy R /V R 
= Kffly . MH yRAW / V R 

- Kffly= cHy GC . V R / MH Yraw 

On tient compte du retard de 6 minutes sur la mesure de cHy GC en 
25 introduisant deux LAGs en sene de 3 minutes chacun sur la valeur du 
modele, et on obtient la formule finale : 

=» Kffly = cHy GC . V R / LAG ( MH yRAW , 3/60, 3/60 ) 

D'une facon similaire a Kffly, on calcule le parametre de correction 
pour la puretd du butene "KfBt" : 
30 Kffly = cBtQ C . V R / LAG ( MBt^ , 3/60, 3/60 ) 
L'algorithme de regulation a pour entrees : 

- les consignes de concentration calculees par Talgorithme maitre; plus 
precisement les consignes pour les rapports de concentration cHy R /cEt R 
"HyEt Sp " et cBt R /cEt R "BtEt sp " (en kg/kg) 

- la consigne pour 1'aUure de marche du precede- FpEt sp , fixee par 
rop^rateur 
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- la consigne pour la concentration en Ethylene cB sp , fixee par 
I'operateur, 

- Ies concentrations calculees par le modele. 

n calcule les consignes pour les debits d'alimentation des reactifs 
FEtgp, FCAgp, FHy sp et FBt sp . Differents algorithmes peuvent etre 
utilises, parmi lesquels le MBPC (Model Based Predictive Control). Us 
peuvent en general se decomposer en un feed-forward base" sur le modele 
inverse, et un feed-back proportionnel a 1'ecart entre le modele direct et les 
consignes de concentration. 
Regulation de 1'alimftntati on en ahvlfen* : 

Feed-forward : valeur pour maintenir la concentration actuelle, basee sur 
Pinversion de la valeur stationnaire de I'equation suivante : 

MEt R = LAG(FEt IN .T Et ,r Et ) ' (kg) 
=» FEtpp = MEtR/ra 

Feed-back : proportionnel a 1'ecart entre la consigne cEtgp, et le modele 

FEtpg m 5 . (cEtgp . V R - MEt R ) 
Consigne : 

FEtgp - FEtpp + FEtpB 
Regulation de Palimenta tion en catalvseiir : 

Feed-forward : valeur pour maintenir la concentration actuelle, basee 
sur I'inversion de la valeur stationnaire de I'equation suivante : 

MCA R = LAG (FCA^ . r CA , r CA ) (kg) 

- FCApp = MCA R /r CA 

Feed-back : proportionnel a I'ecart entre la consigne FpEt sp , et le modele, 
selon la formule suivante : 

FpEt = A^ . MEt R . MCA R / V R 

- FCApg = 5 . (FpEt sp / (A a . MEt R / V R ) - MCA R ) 
Consigne : 

FCA sp = FCApp + FCApb 
Regulation de 1 'alimentation en hvdrogene : 

Feed-forward : valeur pour maintenir la concentration actuelle, basee 
sur 1' inversion de la valeur stationnaire de I'equation suivante : 

MHy R = LAG (FH yiN . r Hy , r Hy ) (kg) 

- FHypp - MH yRAW / r Hy 

Feed-back : proportionnel a 1'ecart entre la consigne de rapport HyEt sp , et 
le modele. 
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FHypg = 5 . (HyEt SP . MEt R - MHy R ) 
Consigne : 

FHy sp = FHypj, + FHypB 
Regulation de 1'alimen tation en hutene (similaire a 1'hydrogene) : 

5 ratEF-MBtRAW^ 

FBtpg = 5 . (BtEt sp . MEt R - MBt R ) 
FBt sp = FBtpp + FBtpg 

Les equations qui precedent resument les equations de 1'algorithme 
esclave. Hies sont executees toutes les 10 secondes par le systeme xaim6- 
10 rique de contiole et de commande (SNCQ. 
L'algorithme ma?tf» 

La Figure 5 illustre le principe de ralgorithme maitre : 

1. son oigane de prediction (baa$ sur un modele direct) pretiit les 
proprietes principales du polymere (MI et MVS); il utilise pour cela la 

15 mesure de la temperature de polymerisation, les predictions pour les 

concentrations dans le reacteur, foumies par le modele esclave, et les 
temps de sejour du PE dans les differents appareils ; 

2. son oigane d'adaptation compare les mesures de MI et de MVS 
effectuees (a la sortie du secheur dgalement), soh par le laboratoire de 

20 mesure toute les 2 h, soit par un analyseur en continu, aux valeurs 

pnSdites par le modele direct, de facon a determiner les 2 parametres 
d'adaptation, qui sont des parametres correctifs, multipUcatif pour le MI 
et additif pour la MVS ; 

3. son organe de controle (base sur un modele inverse) calcule, a paitir des 
25 consignes de MI et de MVS fournie par 1'operateur, les consignes pour 

les concentrations dans le reacteur (rapports Hy:Et et Bt:Et). Comme 
pour ralgorithme esclave, ce calcul est compose d'un feed-forward base 
sur le modele direct, et d'un feed-back proportionnel a I'ecart entre le 
modele direct et les consignes de 1'operateur. 
30 Pour un catalyseur donnd, les proprietes des resines en regime stable 

sont des fonctions de la temperature de polymerisation et des concentrations 
des reactifs. Parmi les differentes equations statiques renseignees dans la 
litterature, on a choisi les equations suivantes : 
log(MI) = hq + aj . T° + a 2 . log (Hy/Et) *+ a 3 . Bt^t 
35 MVS - b 0 + b! . T + bj . (Bt/Et) b3 + b 4 . log(MI) 

Les parametres hq a a 3 et b Q a b 4 sont obtenus par identification en 
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regime stable, pour plusieurs lysines fabriquees avec le meme catalyseur. 

Par ailleurs, les differents appareils rencontres par le polyethylene 
jusqu'au moment ou ses propriety sont mesur&s (reacteur, stripper, 
centrifugeuse puis secheur), peuvent tous etre assimiles en premiere appro- 
5 ximation a des reacteurs paifaitement melanges. 

L'algorithme maitre utilise comme entries les mesures suivantes : 
T R = temperature au reacteur (°Q 
Vstp = volume liquide dans le stripper (obtenu par mesure de niveau) (m^) 
^MES = mesure de MI (melt-index) 
10 MONIES = mesure de MVS (Masse Volumique Sp&ifique) 
. ainsi que les calculs suivants realises par ralgdrithme esclave : 

FpEt = production instantanee de polymere (allure de marche) (kg/h) 



FPEquy = d6bit de PE sonant du reacteur * (kg/h) 

MPE R = masse de PE dans le reacteur (kg) 

15 HyEt R = rapport Hy sur Et dans le reacteur (kg/kg) 

BtEt R = rapport Bt sur Et dans le reacteur (kg/kg) 



Les valeurs instantanees brutes (avant adaptation) de la MVS et du 
logarithme du MI ( n lMI n ) sont calcul6es par : 
MVS^ = b Q + bj . T R + b 2 . (BtEt R ) b3 + b 4 . lML ms 
20 IMIq^s - aQ + a| . T R + a 2 . log (HyEt R ) + a 3 . BtEt R 

En utilisant le theorfeme 2, on calcule les propria moyennes brutes k 
la sortie du reacteur : 
IMIr = LAG (IMIjns, MPE R / FpEt) 
MVSr = LAG (MVS ms , MPE R / FpEt) 
25 En effet : 

- les propriety 1MI et MVS repondent assez bien a une loi de melanges 
lineaire 

- le reacteur boucle peut etre assimite k un reacteur parfaitement melange 

- le debit massique de PE "entrant" (apparaissant) dans le r&cteur est 
30 effectivement FpEt, la quantity de PE qui polymerise k chaque instant 

(allure de marche). 

Proprietes brutes k la mesure : Sachant qu'on a dans le stripper 
environ 500 kg de PE par m^, et si Ton prend rhypothfcse que le stripper 
est un reacteur parfaitement melange, on calcule les proprietes brutes en 
35 sortie du stripper comme suit : 

IMIstp = LAG (IMIr, 500 . Vstp / FPEqut) 
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MVSstp = LAG (MVSr, 500 . Vstp / FPE QUT ) 

Le temps de sejour dans la centrifiigeuse &ant tres court, on peut 
Tignorer. 

Le s&heur est k lit fluide; il contient en permanence environ 
5 1 400 kg de PE. On peut faine l'hypothfcse que le niveau dans le stripper 

change peu, et que le d<5bit qui en sort dgale celui qui y rentre. Des lors, le 

d<3>it de PE entrant dans le secheur vaut FPEquj. On a alors en sortie du 

s&heur, a 1'endroit ou I'6chantiIlon est preleve pour la mesure de propriete, 

les valeurs brutes suivantes : 
10 IMIsh = LAG(lMIstp, 1400/FPEqut) 

MVSsh = LAG (MVSstp, 1400 / FPEq^j) 

Les proprietes aprfcs adaptation sont obtenues en faisant intervenir les 

parametres d'adaptation kMI (parametre multiplicatif) et kMVS (parametre 

additif); les proprtet&i apres adaptation sortie reacteur, stripper et sdcheur 
15 valent done : 

MIr c - kMI . lO 1 ^ 

MVSr c « kMVS + MVSr 

MIstp c = kMI . lolMIstp 

MVSstp^ = kMVS + MVSstp 
20 MIshc = kMI . lolMIsh 

MVSsh^ = kMVS + MVSsh 

Adaptation de I'algorithme maitrp : 

Les mesures de propriety prennem un certain temps k etre realis&s 

(± 5 rain, si analyseur en ligne, ± 1 h si effectuees par le laboratoire). 
25 Pour pouvoir calculer les parametres d'adaptation, il faut done 

resynchioniser (recaler dans le temps) les predictions brutes du modfcle avec 

les mesures. Cela peut par exemple se faire au moyen d'un registre i 

d&alage (appele ici "fonction DELAY") : 

1MI DEL - DELAY (IMIsh, r m ) 
30 MVS DEL = DELAY (MVSsh, tmvs) 

avec : tj^j et rj^yg = ± 5 min. ou ± 1 h, suivant que ia mesure est 

realisfe par un analyseur continu ou par le laboratoire. 

A chaque nouvelle mesure de MI ou de MVS, on recalcule le 

paramfctre d'adaptation brut kMI' ou kMVS', en comparant la valeur 
35 modfele brute resynchronisee a la valeur mesuree : 
" tog ^ (MImes) - 1MI DEL 
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kMVSV = MVSj^es - MVS DEL 

Ces valeurs brutes sont filtr&s afin d'attdnuer les reactions 
precipittes que les perturbations (bruit) €ventuelles des mesures infligeraient 
au precede : 
5 kMI = LAG ( kMT, ± 1 h) 

kMVS = LAG ( kMVS\ ± 1 h) 
Organe de controle : 

L'organe de condyle a pour consignes les valeur MI S p et MVSgp, 
entries par 1'opgrateur. II calcule les consignes pour les rapports des 
10 concentrations dans le reacteur HyEtgp et BtEtgp n&essaire pour obtenir 
rapiderrient les proprietes voulues MIgp et MVSgp. Ce calcul se fait en 2 



1. l'organe de controle calcule, a partir des consignes MI sp et MVSgp 
fournies par l'opdrateur, et des valeurs aprfes adaptation du MI et de la 

IS MVS dans les difterents appareils, les consignes MIi sp et MVSi sp pour 

la production instantanee. Ces consignes instantanees sont composees 
d'un feed-forward, et d'un feed-back proportionnel & T&art entre le 
module direct et les consignes de l'operateur. 

2. les consignes pour les rapports de concentration HyEt sp et BtEt SP sont 
20 alors calculges en inversant I'equation statique utilisee plus haut pour le 

calcul de la valeur instantanee du MI et de la MVS. 

Consignes pour les proprietes instantanees : On compare les 
proprietes sortie secheur aux consignes de proprietes, pour determiner les 
consignes voulues pour les proprietes a la sortie du stripper (on neglige la 
25 centrifiigeuse) : 

MIstp S p = 10 ( log(MISP) + 0.1 . (log(MISP) - log(MIshC))) 

v MVSstpgp = MVSgp + 0.1 . (MVSgp - MVSsh c ) 

De meme, a partir de l'ecart entre ces consignes sortie stripper et 
des valeurs calibrtes au stripper, on calcule les consignes voulues en sortie 
30 du reacteur : 

Mh sp = 10 < log(MIstpSP) + 0.5 . (log(MIstpSP) - iog(MIstpC))) 

MVSr sp = MVSstp sp + 0.5 . (MVSstp SP - MVSstp c ) 

Enfin, & partir de Tecart entre ces consignes sortie reacteur et des 
valeurs calibrees correspondantes, on calcule les consignes voulues pour la 
35 production instantanee : 

MESP = 10 ( log(MKP) + 2 . (log(MIrSP) - log(MIrC))) 
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MYSi sp = MVSr sp + 2 . (MVSr sp - MVSr c ) 

Consignes pour les rapports de concentrations : Les consignes pour 
les rapports de concentration HyEt Sp et BtEt sp s'obtiennent en inversant 
l'equation statique utilisee plus haut pour Ie calcul de la valeur instantanee 
du MI et de la MVS, en substituant aux termes MI et MVS les consignes 
voulues pour la production instantanee et en appliquant le parametre 
d 'adaptation. 

Paitant de : 

log (MBgp- / kMI) = ao + aj . T R + a 2 . log (HyEt sp ) + a 3 . BtEt R 
MVSisp * kMVS = bo + b, . T R + bj . <BtEt Sp )b3 + ^ . 
on obtient : 



a 2 . log (HyEt sp ) = log <MIi sp / kMI) - (ag + a, . T R + a 3 . BtEtp) 
ce qui donne : 

HyEtsp = 10 ( ( ,0g < ^ sp 1 kMI) - aO - al . TR - a3 . BtEtR ) / a2 ) 
15 et: 

02 . (BtEt Sp )W = MVSi Sp -kMVS - ( b Q + b, . T R + b 4 . 1MW, ) 
ce qui donne: 4 m * 

BtEt sp = ( ( MVSi sp - kMVS - bQ - b x . T R - b 4 . IMI^ ) / bj ) W 

Les equations qui precedent resument les equations de ralgorithme 
20 maitre. Elles sont executees toutes les 30 secondes par le SNCC. 

Avec ce precede, U est possible de piloter la polymerisation avec une 
grande precision. En particulier : 

- les proprietes controlees (MI et MVS) sont maintenues au plus pres des 
valeurs voulues, avec une dispersion minime 

25 - les changements de qualite (et done des proprietes MI et MVS) sont 
realises avec rapidite et precision 

- les demanages et arrets de la polymerisation, ainsi que les changements 
d'allure de marche du procede, sont realises d'une facon acceieree, tout 
en maintenant le MI et la MVS fort proches des valeurs voulues. ' 

30 Bien que la methode de regulation selon l'invention ait ete presentee 

a l-aide d'un procede de synthese de polyethylene par polymerisation en 
continu d'ethylene, il est entendu que cette methode de regulation sera de 
fa§on generate efficace pour d'autres precedes de synthese, et en particulier 
pour les precedes presentant 1'une ou plusieurs des caracteristiques 

35 suivantes : 

- une regulation multivariable est necessaire parce que plusieuns variables 
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influent sur r ensemble des propriety k regler; 

- la dynamique du processus est lente : melanges en s<Srie. temps morts 
importants; 

- les mesures des propria sont &hantillonn6es avec une faible frequence, 
5 et/ou bruit&s; 

- la regulation doit etre dynamique, c'est-i-dire valable quelle que soit 
Failure de marche du precede*, ainsi que durant les transitions d 'allure de 
marche et de quality (propriety) du produit k synthetiser; 

- il est interessant d'estimer certaines variables non directement mesurees. 
10 Pour pouvoir etre facilement mise en oeuvre avec les techniques 

pr&entees, il suffit que : 

- les equations statiques du processus soient connues (souvent, elles le sont 
au moins dans une certaine mesure, sinon le proc&te ne pourrait fetre 
maitrise); 

15 - la dynamique du processus puisse etre approchfe par des melanges 
parfaits et des temps morts; 
" les mesures necessaires soient disponibles et d'une quality suffisante (en 
particulier, les debits des nSactifs et les debits traversant les reservoirs 
concerns). 

20 L'utilisation paiticuliere de la fonction LAG d&rite ci-dessus, 

notamment dans les theorfemes 1 et 2, peut naturellement etre etendue k des 
m&hodes de regulation basfes sur une structure autre que celle qui a 6t6 
exposee, qui comprend un algorithme maitre et un algorithme esclave 
distincts. EUe peut par exemple 6tre appliquee dans une methode de 

25 regulation ne comportant qu'un algorithme unique. 
Exemples : 

8 essais de synthese de poly&hylfene (PE) de 4 types differents 
(d6finis par leur MI, MVS, ...) ont 6t6 realises, respectivement en utilisant 
une methode de regulation classique et en utilisant la m&hode de 
30 Tinvention. Le tableau ci-dessous resume les constatations qui ont etd faites 
sur la base de nombreuses mesures de 1'indice de fluidite en fondu (melt 
index) des 8 polymferes obtehus. Cpk designe l'indice de capability centnJ 
du precede. 
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Type 
dePE 


Regulation classique 


Regulation selon 1'invention 


Ecart-type 


Gpk 


Ecart-type 


Cpk 


I 


0.127 


0.909 


0.059 


2.202 


2 




Q.61 




2.0 


3 




0.48 




1.88 


4 




0.64 




1.09 



On constate que l'indice de capabilite Cpk est plus que double grace 
a 1'utUisation de la methode de rinvention, ce qui indique que les proprietis 
sont environ deux fois rnoins dispersees et/ou mieux centres par rapport 
aux consignes imposees. 
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Tableau des abreviations utilises 





a I 


parametres pour l'equation statique du MI (I = 0. a 3) 






*l ' 


parametres pour l'equation statique de la MVS (I = 0 a 4) 






A Et 


activite catalytique pour Methylene (m 3 kg" 1 .!! -1 ) 


5 


CX GC 


concentration de "x" obtenue par la mesure de 1'analyseur 


(kg/m 3 ) 




cx R 


concentration de "x" dans le reacteur 


(kg/m 3 ) 




CX SP 


consigne pour la concentration de "x" dans le reacteur 


(kg/m 3 ) 




Fpx 


defcit massique de polymerisation de "x" (allure de raarche) 


(kg/h) 




FVOUT 


debit volumique sonant du reacteur 


(m 3 /h) 


10 


Fx IN 


debit massique de "x" entrant 


(kg/h) 




Fx OUT 


debit massique de "x" sortant 


(kg/h) 




kd 


constante de desactivation du catalyseur 


(1/h) 




KfBt 


parametre correctif (d'adaptation) pour le butene 






Kffly 


parametre correctif (d'adaptation) pour Thydrogene 




15 


kMI 


parametre correctif (d'adaptation) pour le MI 






kMVS 


parametre correctif (d'adaptation) pour la MVS 






LAG(.) 


fonction filtre passe-bas du ler ordre 








mesure du MI (Melt Index) (indice de fluidite en fondu) 






MIy 


MI (Melt Index) bnit (non calibre) dans "y" 




20 


MIy c 


MI calibre (avec adaptation) dans "y" 






MIy S p 


consigne de MI (calibree) pour "y" 






MVSmes 


mesure de la MVS (Masse Volumique Standard) 






MVSy 


MVS (Masse Volumique Standard) bmte (non calibre) dans 






MVSy c 


MVS calibree (avec adaptation) dans "y" 




25 


MVSy sp 


consigne de MVS (calibree) pour "y" 






Mx RAW 


masse de "x" brute (non calibre) dans le reacteur 


(kg) 




MXy 


masse de "x" calibree (avec adaptation) dans "y" 


(kg) 
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Rx reactivite de X dans le reacteur 

Vy Volume de "y" 

r R temps de sejour dans le reacteur ^ 

. r x temps de sejour pour M x tt dans le nSacteur 

"x" peut representer les constituants suivants : 



Bt 


Butene 


CA 


Catalyseur 


Et 


Ethylene 


Hy 


Hydrogene 


10 Sv 


Solvant 



"y" peut representer les appareils suivants 
r R&cteur de polymerisation 

stp Stripper 
sh Sdcheur $ lit fluide) 
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Legende des figures 

10 reacteur de polymerisation 

1 1 alimentation en reactifs (matiferes premieres), catalyseur, solvant 

12 circuit de refroidissement 
5 14 conduite de soutirage 

16 stripper 
18 condenseur 

20 chromatographe en phase gazeuse 
22 centrifugeuse 
10 24 secheur a lit fluide 

25 solvant et reactifs a recycler 

26 polyethylene 

27 reactifs k recycler 

28 solvant a recycler 

15 30 consignes de propria du polymfere 

31 consigne d'allure de marche du precede 

32 algorithme maitre 

33 algorithme esclave 

34 consignes de concentration 
20 35 consignes de debits entrants 

36 regulation (PID) des debits 

37 mesures 

38 regulation de temperature 

39 feed-forward de temperature 

25 40 dynamique de la polymerisation : cin&ique chimique et bilan matifcre 

41 simulations de Failure de marche et des rapports Hy/Et et Bt/Et 

42 mesures de temperature, de debits et de concentrations 

43 precede regule 

44 analyse d'un dchantillon du produit synthetisfi par le precede 
30 51 mesure de grandeurs liees au d£roulement du precede 

52 mesure de proprietes du polymere 

53 modele direct ; prediction des proprietes a la mesure 

54 comparaison : correction du modfele (adaptation) 

55 algorithme de regulation, base sur le modele inverse (feed-forward + 
35 feed-back) 
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56 consignes pour les grandeurs liees au procedd 

57 mesures et consignes de debits entrants 

58 mesure de la temperature et prediction des concentrations 

59 raodele direct : equations des propriety en fonction des concentitations 

60 prediction des proprietes 

61 consignes pour ies concentrations dans le nfccteur 

62 mesure des debits entrants 

63 mesure des concentrations dans le reacteur 

64 modele direct ; pnkiiction des concentrations sur base du bilan roatiere 

65 predictions de concentration 

66 comparaison ; calcul des parametres d'adaptation 
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RBV-ENDIC ATIQNS 

1 - M&hode de regulation d'un proc&te de synthese d'au moins un 
produit chimique dans un appareillage coxnprenant au moins un nJacteur (R) 
pouvant &re assimite k un reacteur parfaitement melange, dans laquelie une 
5 . ou plusieurs grandeurs de commande (GC) peimettent d'agir sur le 

ddroulement du procedd en vue qu'une ou plusieurs grandeurs lifes aux 
propriety du produit et/ou au deroulement du procede, appel&s grandeurs 
regI6es (GR), soient 6gales a des consignes corespondantes (C GR ), ladite 
m&hode comprenant les Stapes suivantes : 
10 (a) entnfc de consignes concernant les grandeurs rtglees (C QR ) ; 

(b) calcul, au moyen d'un organe de prediction (OP), de predictions des 
grandeurs regimes (P GR ), sur la base de mesures des grandeurs de 
commande du proc6de (M GC ) ; 

(c) utilisation d'un organe de controle (OC) pour calculer des consignes 
15 des grandeurs de commande du proc&te (C GC ), sur base des 

consignes (C GR ) et des predictions (P GR ) des grandeurs reglees ; 

(d) transmission des consignes des grandeurs de commande du procedS 
(C GC ) a des actionneurs, ou a des oiganes de reglage controlant des 
actionneurs, afin d'agir sur le d6roulement du proc&te ; 

20 dans laquelie 1'oigane de prediction (OP) est base sur un modele 

mathematique du proc&Ie, appele modfele direct (M), et est con?u de telle 
maniere qu'on pr&iit la masse d'au moins un constituant (X) dans le 
reacteur (R) par l'6quation : 

M XR " ^ (FxRin • *X • T X ) 

25 dans laquelie : 

- F XRin est le ma ssique du constituant X entrant dans le reacteur R : 

- r x est le temps de sejour de x dans le reacteur, qui vaut 

r x = Mxr / (E Fxdis) 

dans lequel : 

30 " M XR d ^signe la derniere valeur calculee de la masse du constituant 

X pr^sente dans le reacteur R ; 

E Fxdis d6signe la somme de tous les debits massiques Fxdis avec 
lesquels le constituant X disparait du reacteur R ? notamment par 
reaction et/ou par sortie du r&cteur ; 
35 - la fonction y = LAG (u , 7) est la solution de r equation differentielle 
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U = t • ^ ♦ y 

dt 



calculee avec la valeur instantanee de u et de ,, ainsi qu'avec la deraiere 
valeur de y calculee. 

2 - Methode de regulation selon la revendication 1 , dans Iaquelle la 
consigne d'au moins une grandeur regies (C QR ) est corrigee sur la base de 
1'ecart entre la mesure (M GR ) et la prediction (P QR ) de cette grandeur 
reglee, de facon a ce que la regulation soil efficace meme en presence d'une 
erreur dans la prediction de cette grandeur reglee (P GR ). 

3 - Methode de regulation selon la revendication 1, dans Iaquelle le 
modele (M) du precede est pdriodiquement adapte sur la base de Tecart 
entre les predictions (Pq R ) et les mesures (Mq R ) des grandeurs reglees, de 
facon a ce que le modele du precede foumisse des predictions des grandeurs 
re-glees (P GR ) aussi proches que possible des mesures de ces grandeurs 
(Mgr)- 

4 - M^thode de regulation selon la revendication 3, dans Iaquelle les 
mesures (M GR ) des grandeurs reglees interviennent uniquement dans 
I'adaptation eventueUe du modele du precede, et n'interviennent pas 
directement dans le calcul des consignes des grandeurs de commande du 
precede (Cqc). 

5 - Methode de regulation selon Tune des revendications 1 a 4 
appliquee a un precede de polymerisation, comprenant 1'une ou plusieurs 
des etapes supplementaires suivantes : 

- calcul d'une consigne de temperature dans le reacteur en fonction d'une 
ou plusieurs consignes de proprietes du produit ; et transmission de cette 
consigne de temperature a un ou plusieurs actionneurs permettant de 

25 modifier la temperature dans le reacteur ; 

- calcul d'un bilan thermique pour le reacteur, sur base notamment de 
mesures de temperature; utilisation de ce bilan thermique de facon i 
determiner la ^uantite de polymere synthetisg par unite de temps et/ou la- 
productivity du catalyseur et/ou la concentration d'au moins un reactif 

30 dans le reacteur ; 

- calcul de la quantite de chaleur produite par la polymerisation, grace a un 
calcul de la quantite du ou des reactifs qui polymerisent ; determination 



par ce biais de la quantity la chaieur qu'il faut ajouter ou evacuer pour 
maintenir la temperature du reacteur ; utilisation du resultat dudit calcul 
pour ameliorer la regulation de temperature, de fa?on a respecter au 
mieux la consigne de temperature, notamment en cas de modifications de 
. 1'allure de marche. 

6 - M&hode de regulation selon Tune des revendications 
prdcedentes, dans laquelle la propriete Px R d'un constituant "x" dans le 
reacteur R, assimile k un reacteur parfaitement melange, est calculfe comme 
suit : 

Px R = LAG (Px m , Mx R / F XlN ) 
ou "Px" est une propriete d'un constituant "x w , r^pondant de maniere 
substantielle a la loi de melange lineaire Px 1+ ^ = wj . Pxj + w 2 . Px 2 , 
wj et w 2 etant les proportions massiques de deux fractions 1 et 2 de 
propriete Pxj et Px 2 que Ton melange ; 

Px| +2 est la propria de x k sa sortie du reacteur, aprfes melange ; 
Pxjjj est la propriete du constituant "x" k son entree dans le rfacteur R ; 
Mx R est la masse du constituant x dans le reacteur R ; 
Fx IN ^ le roassique du constituant x entrant dans le reacteur R. 

7 - Methode de regulation selon Tune des revendications pi^cddentes, 
comprenant les etapes suivantes : 

- entrfe de consignes relatives k une ou plusieurs proprietes du produit k 
synthetiser, dans un algorithme maitre ; 

- entree de la consigne d'allure de marche du precede, dans un algorithme 
esclave; 

- calcul de consignes de concentration des constituants dans le reacteur k 
l'aide de l'algorithme maitre, en fonction notamment des consignes et des 
mesures de proprietes du produit ainsi que de mesures ou de predictions 
des concentrations des differents constituants dans le reacteur ; 

- transmission des consignes de concentration calculdes par ralgorithme 
maitre comme grandeurs d'entree a ralgorithme esclave; 

- calcul des consignes de debit des constituants entrant dans le reacteur, a 
l'aide de ralgorithme esclave, en fonction notamment de la consigne 
d'allure de marche du precede, des consignes de concentration, et de 
mesures de debit des constituants entrant dans le reacteur, et 

- transmission des consignes de debit calcuiees k l'aide de ralgorithme 
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esclave a un ou plusieurs actionneurs afin de regler les debits des 
constituants entrant dans le reacteur, 
dans laquelle 1'algorithme maitre et/ou 1'algorithme esclave sont mis en 
oeuvre conformernent a Tune des revendications pr&edentes. 

8 - Methode de regulation selon la revendication 7, caractensee en 
ce que 1'algorithme maitre comprend : 

- un organe de prediction, base sur un modele direct du precede permettant 
de fournir une prediction des propri&es du produit synth<5tis<S en fonction 
de mesures et/ou de predictions des concentrations des constituants ; 

- un organe d'adaptation comparant les predictions de proprietes calculees 
par Porgane de prediction a des valeurs effectivement mesurees sur Ie 
produit synth^ et deduisant de cette comparaison des parametres 
d'adaptation, lesdits parametres d'adaptation intervenant comme 
grandeurs d'entree supplementaires dans ledit organe de prediction de 

15 1'algorithme maitre ; et 

- un organe de contrdle, base sur un modele inverse du proc6de\ pour 
calculer, en fonction des consignes et des predictions de propii&es du 
produit a synthetiser, des consignes de concentrations pour 1'algorithme 
esclave, lesdits parametres d'adaptation intervenant egalement comme 

20 grandeurs d'entree supplementaires dans ledit organe de contrdle. 

9 - Methode de regulation selon 1'une quelconque des revendications 
7 ou 8, dans laquelle 1'algorithme esclave comprend : 

- un organe de prediction, base sur un modele direct du precede permettant 
de fournir une prediction des concentrations d'un ou plusieurs des 

25 constituants sur base d'un bilan de matiere dans le reacteur; 

- un organe d'adaptation, comparant les predictions de concentrations 

calculees par le modele direct a des mesures de concentration et 

deduisant de cette comparaison des parametres d'adaptation, lesdits 

paramfctres d'adaptation intervenant comme grandeurs d'entree 

30 supplementaires dans ledit organe de prediction de 1'algorithme esclave ■' 

et 

- un organe de controle, base sur un modele inverse du procede\ pour 
calculer, en fonction de la consigne d'aUure de marche, des consignes de 

concentration calculees par 1'organe de controle de ralgorithme maitre et 
35 des predictions de concentration calculees par 1'organe de prediction de 
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I'algorithme esclave, les consignes pour les debits entrant dans Ie 
nSacteur, lesdits paramfctres d'adaptation intervenant comme grandeurs 
d'entnSe supplemental dans ledit organe de controle de I'algorithme 
esclave. 

10 - M&hode de regulation selon Tune des revendications 
pn$c£dentes, appliquee a la regulation de la synthase en continu de 
polyethylene par polymerisation d'&hylene dans au moins un nSacteur, les 
reactifs comprenant I'&hylene, I' hydrogene et/ou un comonomfere optionnel, 
la ruction de polymerisation ayant lieu en presence d'un catalyseur et une 
partie du contenu du reacteur dtant preleve en permanence ou par 
intermittence. 

11 - M£thode de regulation selon les revendications 9 et 10, dans 
laquelle r organe d'adaptation de 1'algorithme esclave compare les mesures 
des concentrations en Ethylene (Et), hydrogene (Hy) et/ou comonomfere 
optionnel (Bt) aux valeurs pnSdites par 1'organe de prediction de 
ralgorithme esclave, de fa$on a determiner au moins Tun des paramfetres 
d'adaptation suivants : 

a) T activity specifique du catalyseur pour T&hytene "AEt", en kg/h de 
polyethylene par kg de catalyseur et par kg/m 3 d'ethylene ; 

b) I'erreur de gain sur la mesure de debit d'hydrogene "KfHy" ; 

c) la puret6 de I'alimentation en comonomere "KfBV. 

12 - M&hode de regulation selon Tune des revendications 1 a 9, 
appUquee a la regulation de la synthese en continu de polypropylene par 
polymerisation de propylene dans au moins un reacteur, les reactifs 
comprenant le propylene, 1'hydrogene et/ou un comonomfere optionnel, la 
reaction de polymerisation ayant lieu en presence d'un catalyseur et une 
partie du contenu du reacteur £tant preleve en permanence ou par 
intermittence. 

13 - Methode de regulation selon les revendications 9 et 12, dans 
laquelle 1'organe d'adaptation de ralgorithme esclave compare les mesures 
des concentrations en propylene (Pe), hydrogene (Hy) et/ou comonomfere 
optionnel (Et) aux valeurs predites par 1' organe de prediction de 
ralgorithme esclave, de fa?on a determiner au moins 1'un des paramfetres 
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d'adaptation suivants : 

a) ('activity specifique du catalyseur pour le propylene "APe", en kg/h de 

polypropylene par kg de cataJyseur et par kg/ m 3 de propylene ; 

b) Femur de gain sur la mesure de debit d'hydrogene "KfHy" ; ' 
5 c) la purete de Falimentation en comonomere "Kffit". 

14 - Methode de regulation selon la revendication 8, appliquee a un 
proceed de polymerisation, dans laquelle : 

- on mesure p6riodiquement Findice de fluidity en fondu (MI) et/ou la 

masse volumique standard (MVS) du polymere et/ou sa teneur en 
10 comonomere; 

- Forgane de prediction de Falgorithme maitre calcule des predictions 
brutes de MI et MVS en fonction de la temperature au reacteur des 
concentrations dans le reacteur et des temps de sejour dans les diffdrents 
appareils du circuit de polymerisation ; 

- periodiquement, Forgane d'adaptation de Falgorithme maitre : 

- resynchronise les predictions brutes de MI et MVS, en tenant compte 
du temps ecoule entre la prise des mesures de MI et MVS et 
Fobtention du resultat des mesures, et compare les pactions brutes 
de MI et MVS resynchronisees aux mesures de MI et MVS, 

- calcule un parametre d'adaptation multiplicatif kMI applique a la 
prediction brute du MI pour obtenir une prediction calibree du MI 
et ' 

- calcule un parametre d'adaptation adtftif kMVS applique a la 
prediction brute de MVS pour obtenir une prediction calibree de 

25 MVS. 

15 - Methode de regulation selon Fune des revendications 
preceoentes, appliquee a un precede de polymerisation, dans laquelle on 
evalue une ou plusieurs proprietes du polymere en utilisant une technique 
choisie parmi la spectroscopie infrarouge pioche (NIR), la spectroscopie 
infrimouge a transformee de Fourier (FUR) et ia resonance magnetique 
nucleaire (NMR). 

16 - M6thode de regulation selon Fune des revendications 
precedentes, appliquee a un procede de polymerisation, dans laquelle on 
evalue une ou plusieurs proprietes du polymere en appliquant une relation 
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de correlation predtablie aux resultats de mesures r£alisees par spectroscopic 
infrarouge proche (NIR) a plusieurs longueurs d'onde predetermines en 
fonction de la nature du polymere, choisies entre 0,8 et 2,6 mm. 

17 - Proc6d6 de synthase d'un produit chimique dans un appareillage 
5 • comprenant au moins un reacteur pouvant etre assimil£ k un r^acteur 

parfaitement melange, reguld au moyen de la mdthode de regulation selon 
Tune des revendications prec£dentes. 

18 - Dispositif de regulation d'un proc6d6 de synthase d'un produit 
chimique dans un appareillage de synthase comprenant au moins un r&cteur 

10 pouvant 6tre assimil£ k un r&cteur parfaitement melange, dans lequel la 
m&hode de regulation est conforme k Tune des revendications 1 k 16. 
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(Et. Hy, Bt.) 




Fig. 6 
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